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SUMMARY  v  - 

2 

The  perturbations  of  order  J2  in  the  orbital  elements  of  an  earth 
satellite  are  analysed  by  an  elementary  treatment  that  neglects  terms  of  order 
j|  e  .  The  resulting  expressions  are  combined  into  a  triad  of  cylindrical  polar 
coordinates,  defined  by  a  plane  of  fixed  inclination  and  uniform  rotation  rate, 
since  this  leads  to  very  simple  formulae  for  perturbations  in  coordinates. 

Mean  orbital  elements  are  required  and  ane  introduced  in  a  general  manner 
involving  arbitrary  constants.  The  normal  choice  of  constants  is  such  as  to  make 
both  first-order  and  second-order  perturbation  formulae  as  compact  as  possible 
for  the  cylindrical  coordinates,  the  second-order  results  being  expressible  as 

6r  =  ~  -J2  J2  (R2/3)  f  if  cos  4u  +  2  (26  -  3 1 f )  cos  2u}  , 

6u’  =  -  (R/a)2  f  (f  sin  4u  -  (19  -  20f)  sin  2u} 

an<^  3  2  2 

6c  =  -  jg-  ^2  (R  /a)  f  sin  2i  sin  3u  , 

where  R  is  the  earth’s  equatorial  radius;  a  ,  i  and  u  are  the  satellite's 
(mean/  semi-major  axis,  inclination  and  argument  of  latitude,  and  f  is 
sin^  i  . 


Other  aspects  of  the  J2  'main  problem'  are  considered. 
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1  INTRODUCTION 

1 . 1  Previous  papers  on  the  subject 

The  description  of  the  motion  of  an  artificial  satellite  in  the  axi- 

symmetric  field  due  to  the  lover-order  zonal  harmonics  of  the  geopotential  has 

been  recognized  as  the  'main  problem'  in  the  theory  of  satellite  orbits,  and  many 

solutions  have  been  published.  The  first  RAE  paper'  appeared  in  the  same  month 

2 

as  the  first  satellite  and  was  succeeded  by  the  theory  which  Merson  developed  as 

3 

a  basis  for  the  RAE's  first  proper  orbit  determination  program  .  Merson  took  this 
first  theory  as  far  as  J,  (results  for  an  arbitrary  J.  were  obtained  by  the 

4  °5  z  6 

present  author  )  and  then  compared  it  with  the  theory  of  Kozai  (taken  to  ) 

that  pioneered  the  orbit  analysis  of  the  Smithsonian  astrophysical  Observatory. 

In  a  second  theory^,  Merson  developed  the  Kozai  approach  as  the  basis  for  a  new 

8  9 

orbit  determination  program,  PROP,  that  *  is  still  in  use.  The  paper  by  Kozai 
may  be  regarded  as  one  of  the  two  classic  American  papers  (its  approach  is 

summarized  in  a  recent  text-book'^),  the  other  being  the  paper  of  Brouwer'' 
published  at  the  same  time.  The  Brouwer  approach  is  followed  in  the  book  of 

12 

Brouwer  and  Clemence  ,  whilst  one  of  the  first  text-books  to  give  results  (for 

1 3 

J2  only)  was  that  of  Sterne  . 

The  early  papers  normally  provided  long-term  solutions  of  the  'main 
problem',  ie  gave  formulae  for  perturbations  in  mean  orbital  elements,  to 

n  o 

0(0  ,  where  J.  ■  0(J,)  if  l  >  2  ,  but  only  took  short-periodic  perturbations 
to  0(J2)  ,  neglecting  them  entirely  for  higher-order  harmonics.  It  was  found  by 
Vinti1^-1^  that,  by  use  of  spheroidal  coordinates,  a  complete  solution  could  be 
found  for  a  field  including  J,  ,  but  the  higher-order  even  harmonics  had  to  be 
powers  of  J2  (viz  =  -J2,  -  J2,  etc),  so  the  general  problem  remained. 

2 

The  first  paper  to  take  short-periodic  perturbations  to  0(J»)  seems  to 

18  ^  19 

have  been  that  of  Petty  and  Breakwell  .  Kozai,  in  a  comprehensive  development 

2  .  3 

to  0(J2)  for  short-periodic  perturbations  and  0(J2)  for  long-term  perturba¬ 
tions,  changed  his  approach  to  the  method  of  von  Zeipel  that  had  been  used  by 
Brouwer.  Aksnes  followed  Hori  in  basing  his  approach  on  the  use  of  Lie 
series;  by  reference  to  a  suitable  intermediate  orbit,  he  thereby  obtained 
results  that  were  equivalent  to  Kozai 's  but  much  more  compact.  Much  more  elemen¬ 
tary  (and  therefore  comprehensible)  methods,  just  based  on  the  planetary  equations 

2 

of  Lagrange,  were  adopted  by  the  French,  the  short-periodic  perturbations 

22  3  ^  23 

being  given  by  Bretagnon  and  the  J2  long-term  effects  by  Berger 

Solution  of  the  problem  to  get  a  further  order  -  short-periodic 

perturbations  to  0(J2)  and  long-term  effects  to  0(J2)  -  was  achieved  by 
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Deprit  and  Rom  by  the  application  of  computer  algebra  to  a  method  based  on  Lie 

series.  Their  results  are  only  valid  for  moderately  small  values  of  the  orbital 

eccentricity,  e  ,  since  they  are  expressed  as  truncated  power  series  in  e  . 

Techniques  based  on  computer  algebra  are  obviously  very  powerful,  and  two  papers 

have  recently  been  published  that  give  main-problem  solutions  involving  vast 

25 

numbers  of  terms.  Berger  and  Walch  cover  harmonics  as  far  as  ;  they  only 

go  to  second  order  in  short-periodic  perturbations  and  (basically)  third  order  in 

long-term  effects,  but  this  means  the  inclusion  of  secular  terms  in  , 

for  example,  that  originate  as  fourth-order  couplings  between  JL  and  J_  ,  so 

26  ^ 

that  a  large  number  of  combinations  are  involved.  Kinoshita  covers  only  the 

classical  main-term  harmonics,  te  to  only,  but  goes  to  third  order  for  short- 

periodic  perturbations  and  to  fourth  order  for  long-term  effects;  the  perturba- 
21 

tion  method  of  Hori  is  adopted  and  the  resulting  accuracy,  for  satellite 
motion  in  a  low-eccentricity  orbit,  is  better  than  1  cm  over  a  month. 

Much  of  the  literature  cited  is  extremely  difficult  to  follow,  even  for  a 
reader  with  considerable  experience  of  orbital  analysis,  mainly  because  of  the 
extreme  sophistication  of  the  methods  used,  but  probably  also  because  of 
unfamiliarity  with  the  notation.  Without  understanding  the  methods  it  is  diffi¬ 
cult  to  interpret  the  results  -  to  pick  out  (say)  a  dominant  set  of  terms  for 
2 

(K^)  perturbations  for  a  near-circular  orbit  becomes  virtually  impossible. 

This  is  particularly  true  if  expansions  are  given  in  terms  of  a  quantity  q  , 
defined  as  /l  -  e^  (this  quantity  is  denoted  by  q  in  the  present  paper). 

1 .2  The  present  paper 

In  many  applications,  involving  satellites  in  orbits  of  low  eccentricity, 

.  2 

there  are  significant  perturbations  of  order  J~  ,  including  terms  of  short 

.  2  .  .f 

period,  but  terms  m  ^  e  are  entirely  negligible.  The  significant  terms  may 
be  obtained  quite  easily,  by  an  elementary  iteration  on  Lagrange's  planetary 
equations,  if  a  first-order  solution  to  e)  is  already  available.  Further¬ 

more,  when  an  appropriate  set  of  non-singular  elements  is  used,  the  results  can 
be  expressed  very  compactly.  This  is  the  approach  of  the  present  paper,  the 
particular  set  of  non-singular  elements  being  (as  defined  in  section  2.1) 
a,  i,  n,  5,  p  and  U  . 

The  paper  makes  two  contributions  to  the  solution  of  the  'main  problem', 
and  discusses  various  aspects  of  the  philosophy  of  perturbation  theory.  The 
first  contribution  concerns  the  derivation  of  formulae  for  perturbations  in 
position,  not  merely  orbital  elements.  If  the  satellite's  position  can  be 

j 
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expressed  by  three  suitable  coordinates,  then  the  number  of  perturbation  formulae 
required  is  immediately  halved.  (In  certain  applications,  for  example  the 
processing  of  Doppler  data,  formulae  for  three  velocity  components  will  also  be 
required,  but  these  can  be  obtained  from  the  position  formulae,  essentially  by 
differentiation.)  The  reduction  will  be  all  the  greater  if  the  formulae  for 
position  are  simpler  than  those  for  certain  of  the  elements,  as  is  clear  from 
Kozai's  original  paper^.  By  introducing  r  (geocentric  distance)  and  u  (argu¬ 
ment  of  latitude),  Kozai  was  able  to  reduce  long  expressions  for  the  (first-order) 
short-periodic  perturbations  in  a,  e,  o>  and  M  (defined  here  in  section  2.1) 
into  compact  expressions  for  the  perturbations  in  r  and  u  ;  the  expressions 
for  the  perturbations  in  i  and  ft  are  much  shorter  and  he  left  these  alone. 

In  developing  the  Kozai  approach  for  use  with  PROP,  Merson  presented  six 

perturbation  formulae  that  are  extremely  compact  and  totally  free  of  singularity, 

.  9 

and  the  present  author  modified  them  slightly  to  reduce  non-linear  effects. 

There  is  a  different  way  to  develop  the  Kozai  approach,  however,  and  in  a 
27 

recent  study  of  the  perturbations  of  a  Navstar  satellite  due  to  a  general 
harmonic  I  presented  results  in  what  I  called  the  'RLC  system'  of  mutually 

orthogonal  directions,  where  R  refers  to  r  (as  above),  L  to  an  along-track 
quantity  l  that  only  exists  at  the  differential  level  (and  is  such  that  Si 
is  the  product  of  r  with  one  of  the  six  perturbations  of  Merson,  viz 

6u  +  6ft  cos  i),  and  C  to  a  cross-track  displacement  (such  that  6c  is  the 

.  .  .  .  .  2 
product  of  r  with  6i  sin  u  -  6ft  sin  i  cos  u) .  Ref  27  does  not  cover 

effects,  however,  and  the  RLC  system  is  not  really  suitable  for  effects  beyond 

the  first  order  in  any  J  ,  because  an  integral  l  is  not  defined;  but  it  is  a 

small  step  -  taken  here  in  section  2.5  -  to  a  system  that  suitable.  The 

system  amounts  to  a  set  of  cylindrical  polar  coordinates  defined  in  reference  to 

a  steadily  rotating  mean  orbital  plane,  and  the  formulae  for  second-order  short- 

periodic  perturbations  in  these  coordinates  constitute  the  main  results  of  this 

paper . 

The  paper's  other  contribution  to  the  main  problem  relates  to  the  defini¬ 
tion  of  'mean'  orbital  elements,  and  hence  to  the  'mean  orbital  plane'  just 
referred  to.  The  concept  of  a  mean  element  underlies  all  approaches  to  the  main 
problem,  but  there  is  an  ambiguity  in  every  such  element  in  that  the  short- 
periodic  perturbation  (which,  added  to  the  mean  element,  gives  the  uniquely 
defined  osculating  element)  is  arbitrary  by  any  quantity  that  is  free  of  short- 
periodic  perturbation.  Thus  if  £  is  one  of  the  six  orbital  elements  and 
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(we  later  make  such  expansions  in  terms  of  a  quantity  R  ,  related  to  Jg  > 

rather  than  Jj  itself),  then  is  arbitrary  to  the  extent  of  a  quantity  that 

could  conveniently  be  denoted  by  and  is  itself  a  function  of  the  six  mean 

elements.  Furthermore,  if  j  >  1  ,  each  is  necessarily  a  function  of  the 

k.  for  smaller  j  . 

J? 

Most  papers  on  the  main  problem  do  not  explain  how  their  mean  elements  are 
defined,  nor  indeed  admit  that  there  is  any  problem.  This  makes  it  difficult  to 
compare  the  results  of  different  theories,  and  -  worse  -  hard  to  interpret  the 
published  elements  of  actual  satellites.  The  approach  of  the  present  paper  is  to 
face  up  to  the  problem  by  introducing  the  constant-like  k-quantities  explicitly 
into  formulae,  so  that  general  results  can  be  given  without  a  pre-empted  inter¬ 
pretation  of  particular  mean  elements.  Individual  first-order  k^  (effectively 
the  k^  with  appropriate  'scaling')  are  introduced  in  section  3  (subsequent  to 
further  discussion  of  mean  elements  in  section  2.4),  the  most  important  being 

k  ,  k .  and  k 
a’  i  n 

There  are  three  possible  criteria  which  can  underlie  the  choice  of  the 
k-quantities  and  which  (tacitly  at  least)  are  applied  in  the  literature.  First, 
the  k's  may  be  chosen  so  as  to  make  the  expressions  for  positional  perturba¬ 
tions  as  simple  as  possible;  this  is  the  philosophy  adopted  here.  Secondly,  they 
may  be  chosen  such  that  time  averages  of  the  short-periodic  perturbations  are 
strictly  zero;  this  amounts  to  making  the  k's  zero  if  perturbations  are  expres¬ 
sed  in  terms  of  mean  anomaly  (M)  rather  than  true  anomaly  (v) .  Thirdly,  the  k's 
may  be  set  to  zero  when  the  perturbations  are  expressed  in  terms  of  v  ;  higher- 
order  formulae  will  then  be  simpler  than  with  the  second  choice. 

This  brings  us  to  another  important  point.  Early  formulations  of  the 

short-periodic  perturbations  ( eg  Refs  6  and  13)  were  in  terms  of  v  (or  u  , 

where  u  -  v  =  w  ,  the  argument  of  perigee)  and  hence  expressible  in  closed  form. 

Later  papers,  and  in  particular  those  based  on  computer  algebra,  switched  to  M  , 

for  reasons  that  are  not  entirely  clear,  the  result  being  expressions  that  are 

necessarily  truncations  of  power  series  in  e  .  The  present  paper,  even  though 

2 

it  truncates  before  the  OC^  e)  terms,  uses  v  (or  rather  u  );  there  is  no 
obvious  reason  why,  in  terms  of  v  ,  the  main  problem  should  not  be  solved  to 
further  order  without  a  requirement  for  infinite  series  in  e  (but  see  also  the 
remarks  of  Refs  24  and  28;  the  difficulty  lies  in  evaluating  JvdM). 

Section  9  of  this  paper  is  devoted  to  complete  expressions  for  the  second- 

order  long-term  variation  of  the  orbital  elements,  divided  into  secular  and  long- 

2 

periodic  perturbations  in  the  usual  way.  When  (K^  e)  effects  can  be  neglected. 
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no  long-periodic  perturbations  arise  (the  meaning  of  'second-order'  in  connection 
with  these  perturbations  is  discussed),  but  (as  already  indicated)  it  makes  sense 
to  take  long-term  perturbations  to  a  higher  order  than  short-periodic  perturba¬ 
tions,  unless  the  satellite's  motion  is  only  to  be  represented  for  a  matter  of, 
say,  at  most  a  revolution  or  two.  For  completeness,  also,  the  untruncated 
expressions  for  the  first-order  perturbations  in  positional  coordinates  are 
given  in  section  8,  since  section  3  takes  them  only  to  CK^  e)  . 

The  results  obtained  for  perturbations  in  orbital  elements  are  compared 

22  25  .  26 

with  those  given  by  Bretagnon  ,  Berger  and  Walch  ,  and  Kinoshita  ,  all  dis¬ 
crepancies  being  resolved.  The  formulae  for  second-order  positional  perturba¬ 
tions  have  been  checked  by  comparing  computed  results  with  those  given  by 
numerical  integration. 

2  BACKGROUND 

2.1  Osculating  elements 

The  standard  osculating  elements  of  an  elliptic  orbit  are  a  (semi-major 
axis),  e  (eccentricity)  ,  i  (inclination)  ,  fl  (right  ascension  of  the  node)  , 

<d  (argument  of  perigee)  and  M  (mean  anomaly)  .  Since  e  is  here  assumed  to  be 
small,  we  shall  also  require  the  non-singular  elements  £  ,  n  and  U  ,  defined 
by 

£  =  e  cos  u>  ,  (1) 

n  =  e  sin  w  (2) 

and 

U  =  M  +  u  .  (3) 

The  mean  motion,  n  ,  though  defined  by 

2  3  ... 

n  a  -  v  ,  (4) 

where  v  is  the  earth's  gravitational  constant,  will  often  be  treated  as  an 
element  in  its  own  right.  Finally,  it  is  convenient  to  introduce  f,  h,  p  and 
q  ,  where 

f  ■  sin2  i  ,  (5) 


1  -  lif  , 


p  -  a  ( 1  -  e  ) 

q  -  0  -  e V  . 
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2.2  Assumed  potential 

As  we  are  only  concerned  with  the  Jj  term  of  the  earth's  disturbing 
potential,  U  ,  the  potential  is  taken  as  p/r  +  U  ,  where 

U  =  '  r  J2  (f)2  P2(ain  e)  ;  (9) 

here  r  is  the  radial  distance  of  the  satellite  from  the  centre  of  the  earth, 

8  is  its  geocentric  latitude,  R  is  the  earth's  equatorial  radius,  and  is 

the  usual  Legendre  polynomial.  Thus  we  are  neglecting  drag,  lunisolar  attraction, 
etc,  as  well  as  the  other  harmonics  of  the  real  geopotential. 

To  eliminate  8  ,  we  introduce  the  argument  of  latitude,  u  ,  of  the 
satellite,  given  by 


U  =  V  +  to  , 

where  v  is  the  true  anomaly  (see  section  2.5);  thus, 

sin  6  =  sin  i  sin  u 


(10) 


and  hence 


U 


=  i  H  J2  3  |2  -  3f  (1  -  cos  2u)J 
r 


(ID 


It  is  easy  to  eliminate  r  ,  since 


—  =  1  +  e  cos  v 

r 


(12) 


2.3  Lagrange's  planetary  equations 

Rates  of  change  of  osculating  elements  may  be  expressed  exactly  by 
Lagrange's  planetary  equations  *12,13^ 


• 

2  ZU 

a 

na  3M 

f 

• 

1 

1  2 

W 

e 

“  2  * 

3M 

na  e 

• 

• 

cosec 

i  1 

2 

cos  1 

na  q 

3gu 


W) 

an;  ’ 


(13) 


(14) 


(15) 
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cosec  i  3 U 

2  IT 


na  q 


(16) 


and 


i  ja  a u_ 

2  e  3e 
na 


cot 

i  W 

q 

3l 

n  +  a  , 


(17) 

(18) 


where 


1 

2 

na 


qf.  W_ 

e  Be 


+ 


(19) 


Equation  (19)  effectively  defines  o 
(t  *  0)  ,  on  the  basis  that 


M  = 


as  the  'modified  mean  anomaly  at  epoch' 


(20) 


It  will  be  observed  that  U  ,  where  U  is  given  by  (3) ,  is  not  quite  free  of 
3f//3e  ,  since  3t//3e  is  multiplied  by  different  powers  of  q  in  (17)  and  (19). 
Thus  though  U  is  a  better  element  than  M  for  use  when  e  is  small,  it  can 
complicate  matters,  rather  than  simplify  them,  when  a  full  analysis  is  required. 

The  partial  derivatives  required  in  Lagrange's  equations  can  be  obtained 
from  (11),  bearing  in  mind  that 

3v  m  sin  v  (2  +  e  cos  v) 

3e  2 

q 

and 

2  2 
3v  a  q  (1  +  e  cos  v) 

3M  2  5  • 

r  q 


It  is  also  convenient  to  introduce,  for  general  use  throughout  the  paper,  the 
quantity  K  ,  where 


K 


(21) 


so  that  we  may  re-express  U  as 


U  -  -r  —  (2h  +  3f  cos  2u)  (1  +  e  cos  v)3  .  (22) 

6  P 

(It  has  to  be  remembered,  from  time  to  time,  that  K  is  not  constant;  thus 
3K/3a  -  and  less  importantly  3K/3e  -  is  not  zero.) 
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Then  exact  expressions  for  the  rates  of  change  of  the  elements  are: 


a  = 


-  ^ 4 f  sin  2u  +  e  ^4h  sin  v  -  f  sin  (u  +  <o)  +  5f  sin  (2u  +  v)j|  , 

.  (23) 


Kn 

4q3 


4h  sin  v  +  f  sin  (u  +  to)  +  7f  sin  (2u  +  v) 

+  [4h  sin  2v  +  12f  sin  2u  +  5f  sin  2(u  +  v)  -  f  sin  2io]|  , 

.  (24) 


Kn  sin  2i  /p\  .  „ 

- - -  f  sin2u  , 


2q' 


(25) 


a  = 


Kn  cos 


-(?) 


(1  -  cos  2u) 


(26) 


e  to 


Kn 

4q3 


and 


e  a  =  - 


Kn 

4q: 


|4h  cos  v  -  f  cos  (u  +  to)  +  7f  cos  (2u  +  v) 

+  { e  [4h  cos  2v  -  2  (4  -  7f)  cos  2u  +  5f  cos  2(u  +  v) 

+  f  cos  2a)  +  2(6  -  7f)]j 
.  (27) 

^4h  cos  v  -  f  cos  (u  +  to)  +  7f  cos  (2u  +  v) 

+  je  [4h  cos  2v  -  18f  cos  2u  +  5f  cos  2(u  +  v) 

+  f  cos  2co  -  12hj|  . 


(28) 


The  last  two  equations  are  written  with  e-factors  on  the  left-hand  side  to 
emphasise  their  singularity.  The  sum,  b  +  <o  ,  is  free  of  singularity,  and  will 
be  used  in  section  4.7,  but  as  an  exact  combination  it  does  not  simplify,  because 
of  the  difference  in  powers  of  q  already  remarked. 

•  • 

Exact  expressions  for  £  and  q  are  easy  to  obtain,  since 


and 


£  =  e  cos  u  -  e  (o  sin  w 


q  =  e  sin  to  +  e  to  cos  (o 


(29) 

(30) 


J 


1  uu 
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-  -5^  ^ |  |(4h  +  f)  sin  u  +  7f  sin  3u 

+  £e  [3(4  -  5f )  sin  w  +  (8  -  7f)  sin  (u  +  v) 

-  (4  -  1 3f )  sin  (2u  +  w)  +  5f  sin  (3u  +  v)]| 


*  =,  -SlL  ^  |(4h  -  f)  cos  u  +  7f  cos  3u 

+  [e  [(12-  1 3 f )  cos  oj  -  5f  cos  (u  +  v) 

-  (4  -  1 2f )  cos  (2u  +  u)  +  5f  cos  (3u  +  v)]j  .  (32) 

2.4  Mean  elements 

Let  t,  be  an  osculating  element,  ie  it  stands  for  any  one  of  a,  e,  i,  ft, 
u,  a,  M,  £,  n  and  U  .  Its  variation  is  made  up  of  a  long-term  effect,  not 
directly  related  to  the  orbital  period  ( ie  it  is  independent  of  u  and  v  ),  and 
of  a  series  of  short-periodic  terms  (period  a  submultiple  of  the  orbital  period). 
The  combination  of  the  short-periodi ;  terms  is  the  short-periodic  perturbation, 
(5£  ,  removal  of  which  from  £  gives  a  mean  element,  £  ,  which  for  most  purposes 
(because  its  variation  can  be  plotted  easily  and  accurately  over  long  periods  of 
time)  is  more  useful  than  the  osculating  e  .t;  thus 

S  =  l  +  .  (33) 


The  variation  of  £  in  general  has  two  components,  a  purely  secular  effect, 

which  is  zero  for  several  of  the  elements,  and  a  long-periodic  effect  related  to 

the  secular  variation  of  w  .  Now,  to  the  first  order  in  ^  there  is  no  long- 

periodic  effect  (the  meaning  of  ’order'  in  relation  to  long-periodic  perturbations 

is  discussed  in  section  9),  since  long-periodic  effects  are  induced  by  the  first- 

order  variation  of  u>  and  hence  only  appear  at  higher-than-f irst  order;  at 

2 

second  order,  though  there  is  an  effect,  it  is  e)  ,  which  (until  section  9) 

we  are  neglecting.  Hence  we  may  (normally)  write 

• 

l  =  lQ  +  l  t  ,  (34) 


where  £  vanishes  when  c  is  a,  e  or  i  .  It  will  sometimes  be  useful  to 
write 
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C  -  C0  +  ^  .  (35) 

where,  from  (33)  and  (34) 

AS  =  s  t  +  AS  .  (36) 

We  shall  be  expressing  6s  as  a  polynomial  (truncated  to  a  quadratic)  in 
K  ,  where  K  is  K  ,  as  given  by  (21),  with  p  replaced  by  p  .  It  is  evident 
that  the  'coefficient'  of  each  'term'  of  the  polynomial  is  itself  a  function  of 
the  (mean)  orbital  elements,  possibly  involving  many  terms,  and  that  such  a 
function  is  arbitrary  to  the  extent  of  a  u-independent  term,  ie  that  mean  elements 
are  only  defined  to  within  selected  'constants'.  Choice  of  a  set  of  constants  is 
a  matter  of  convenience,  and  published  papers  do  not  always  make  clear  what  parti¬ 
cular  constants  have  been  chosen,  or  the  motivation  for  the  authors'  choice. 

As  indicated  in  section  1 ,  the  present  paper  generalizes  the  approach  to 
mean  elements  by  explicitly  introducing  arbitrary  'constants'  (k&,  kg  etc) 
into  the  expressions,  Cj  ,  that  occur  as  coefficients  of  K  in  the  various  6c  . 
At  the  second-order  level,  similarly,  we  meet  k2a  etc. 

2.5  Satellite  position 

The  significance  of  a  set  of  mean  elements,  c  »  is  perhaps  best  understood 
by  seeing  how  they  can  be  used  directly  (in  conjunction  with  perturbation 
formulae)  to  generate  a  satellite's  position  coordinates,  viz  x,  y  and  z  . 

The  standard  algorithm  for  x,  y  and  z  ,  given  the  osculating  elements 
a,  e,  i,  n,  w  and  M  ,  is  as  follows: 

(i)  the  eccentric  anomaly,  E  ,  is  found  by  solving  Kepler's  equation 

E  -  e  sin  E  =  M  ;  (37) 


(ii)  the  true  anomaly,  v  ,  is  found  from  one  of  the  two  equivalent 

formulae  (apart  from  an  ambiguity  of  quadrant  in  the  first  formula) 


and 


tan  v 


q  sin  E 
cos  E  -  e 


tan  Jv 


(H i)*  tan  iE 


(38) 

(39) 


100 
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(iii)  u  is  obtained  from  (10),  and  r  from  either  (12)  or  the  equivalent 
formula 


r  =  a  (1  -  e  cos  E) 


(40) 


(iv)  x,  y  and  z  are  obtained  from  the  double  coordinate  transformation 
expressed  by  the  matrix  formula 


(x 


2)  “  Rj  (~i)  (r  cos  u  r  sin  u 


T 

0)  * 


where  R^(0)  describes  rotation  about  the  \th  axis,  so  that 


R,(6)  = 


V 

also  T  ,  in  (41),  denotes  transposition. 
There  is,  of  course,  a  preliminary  step  given  by 


(41) 


M  =  Mq  +  nt  ,  (42) 

if  the  mean  anomaly  is  only  available  at  epoch,  though  appropriate  for  use  at 
time  t  ,  and  it  may  be  necessary  to  recover  e,  u>  and  M  from  non-singular 
elements  £,  n  and  U  defined  by  equations  (1)  to  (3). 

Now  suppose  the  starting  point  is  a  set  of  mean  elements  at  epoch,  with 
secular  rates  and  short-periodic  perturbations  known.  Then  each  £  (osculating) 
can  be  obtained  from  (35)  before  operating  the  algorithm,  but  there  is  an  alter¬ 
native  procedure  which  is  more  satisfactory  since  it  requires  only  four 
6-expressions  instead  of  six.  In  this  procedure,  which  was  given  by  Kozai  , 
quantities  6r  and  6u  ,  representing  short-periodic  perturbations  in  r  and  u  , 
appear  instead  of  6a,  6e,  6u>  and  6M  .  Steps  (i)  to  (iii)  of  the  standard 
algorithm  now  operate  on  the  mean  elements,  leading  in  turn  to  E,  v,  u  and  r  , 
after  which,  as  an  additional  part  of  step  (iii),  r  and  u  are  derived  as 
r  Mr  and  u  +  6u  respectively.  The  operation  of  step  (iv)  is  exactly  as 
before. 

But  Kozai's  modification  of  the  standard  algorithm  for  x,  y  and  z  can  be 
extended,  to  its  logical  conclusion,  so  that  only  three  6-expressions  are  required. 
Two  of  the  perturbed  quantities  are  slightly  modified  forms  of  r  and  u  ,  and 
may  thus  be  conveniently  written  r'  and  u'  ,  whilst  the  third  represents  the 


I 
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displacement  of  the  satellite  from  its  'mean  orbital  plane',  ie  from  a  plane  of 
fixed  inclination,  i  ,  that  rotates  at  a  fixed  rate,  12  .  This  displacement  is 
in  the  'cross-track  direction'  and  may  be  conveniently  denoted  by  c  ,  though  we 
cannot  distinguish  the  perturbation  6c  from  c  itself,  since  the  unperturbed 
value  is  zero.  The  positive  direction  for  c  is  such  that  the  (positive) 
direction  of  orbital  motion  is  given  by  the  right-hand  screw  rule,  the  first- 
order  formula  for  5c  being  in  consequence 

6c  =  r  (6i  sin  u  -  6ft  sin  i  cos  u)  ,  (43) 

as  given  in  Ref  27 . 

This  approach  gives  a  complete  first-order  representation  of  6i  via  6c  , 
but  only  represents  the  6ft  sin  i  component  of  5ft  .  The  other  component  has 
to  be  incorporated  with  6u  ,  the  combined  effect  being  6u'  and  given  to  first 
order  by 

6u'  =  6u  +  6ft  cos  i  .  (44) 

The  difference  between  r'  and  r  is  zero  to  first  order. 

The  operation  of  step  (iv)  is  different  with  this  approach,  since  i  and 
ft  ,  not  i  and  ft  ,  must  be  used  in  the  coordinate  transformation.  Thus  (41)  is 
replaced  by  the  formula 

(x  y  z)T  =  (— Q)  Rj(-i)  (r'  cos  u'  r'  sin  u'  c)T  .  (45) 

It  is  clear  that 

r'2  +  c2  =  r2  ,  (46) 

2  2  2 

since  they  are  equal  to  the  same  thing.  Viz  x  +  y  +  z  .  It  is  from  (46)  that 
we  see  the  first-order  identity  between  r’  and  r  ;  if  6c  can  be  kept  to 

0  o  o  o  / 

0(Ke,  K  )  ,  indeed,  then  r'  -  r  will  be  0(K  e  ,  K  e,  K  )  and  hence  (for  our 

purpose)  entirely  negligible.  In  section  3  we  shall  see  that  5c  can  be  kept  to 

2 

0(Ke,  K  )  by  appropriate  choice  of  k.  and  k^ 

l  u 

The  quantities  r',  u'  and  c  form  a  set  of  cylindrical  polar  coordinates 
relative  to  the  steadily  rotating  mean  orbital  plane.  Thus  the  directions  in 
which  r'  increases,  u'  increases  and  c  increases  are  orthogonal,  correspon¬ 
ding  to  the  axis  directions  for  the  RLC  system  of  coordinates  previously  used  by 
27 

the  author  .  The  trouble  with  that  system  is  that  the  'L'  coordinate  is  only 
defined  differentially  (as  remarked  in  section  1)  on  the  basis  that 
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=  r  6u'  .  (47) 

2.6  Formulae  connecting  mean  and  true  anomaly 

Formulae  connecting  M  and  v  will  be  required  throughout  the  paper. 
Suitable  formulae  are  collected  here  for  reference.  They  are  valid  for  both 
osculating  elements  and  mean  elements. 

The  expression  for  the  difference  between  M  and  v  is  known  as  'the 
equation  of  the  centre'.  The  following  is  an  exact  form  of  the  equation,  expres¬ 
sed  as  a  double  series  of  Bessel  functions,  as  given  by  page  77  of  Ref  12: 


v  =  M  +  2 


t  ‘J7JS|vj*>  ♦  i(rhf)k  *  Wje)]l 


(48) 


J  =  l 


k= 


This  leads  to 


ry  q  y 

v  »  M  +  2e  sin  M  +  1  Je  sin  2M  -  j^re  (3  sin  M  -  13  sin  3M)  +  0(e  )  , 


(49) 


of  which  the  equivalent  v-series  result  is 

ry  o  y 

M  =  v  -  2e  sin  v  +  |e  sin  2v  -  ie  sin  3v  +  0(e  ) 


(50) 


From  (49)  we  can  derive  formulae  for  sin  v  and  cos  v,  viz 

2 

sin  v  *  sin  M  +  e  sin  2M  -  ie  (7  sin  M  -  9  sin  3M) 

-  ie3(7  sin  2M  *  8  sin  4M)  +  O(e^)  (51) 

and 

2 

cos  v  =  cos  M  -  e  (1  -  cos  2M)  -  lie  (cos  M  -  cos  3M) 

-  lie3  (cos  2M  -  cos  4M)  +  0(eS  .  (52) 

Finally,  we  generalize  (51)  and  (52)  to  0(e)  expressions  for  the  sine  and 
cosine  of  (jv  +  ku>)  .  Allowing  T  to  denote  either  one  of  these  functions, 
we  have 

T(jv  +  kw)  -  T(jM  +  k<o)  +  je  T^(j  +  1)M  +  kwj  -  je  T  |(j  -  1)M  +  ku>^  +  0(e2)  . 

.  (53) 

We  can,  of  course,  replace  M  by  v  ,  if  desired,  in  the  two  je-terms.  Also, 
by  taking  j  -  k  ,  we  can  derive  the  relations  between  the  sines  and  cosines  of 
ku  and  kU  . 


16 


3  FIRST-ORDER  RESULTS 

Complete  first-order  perturbation  formulae  may  be  expressed  in  closed  form, 
since  powers  of  e  beyond  the  third  do  not  arise.  To  integrate  equations  (23) 
to  (28),  (31)  and  (32)  with  respect  to  t  ,  we  transform  to  v  as  independent 
variable  by  use  of  the  relation  (for  an  unperturbed  orbit) 


This  leads  to  an  expression  for  dc/dv  ,  for  each  element  £  ,  in  terms  of  a 
power  of  (p/r)  and  trigonometric  terms  in  v  and  w  .  The  power  of  (p/r) 
can  be  expanded  by  use  of  (12),  after  which  the  expression  for  d?/dv  can  be 
integrated  with  respect  to  v  ,  all  elements  being  held  constant. 

The  results  have  been  given,  in  various  equivalent  forms,  by  a  number  of 
authors,  eg  in  Refs  10,  12  and  13.  They  are  given  again,  in  a  compact  form,  here 
for  completeness.  Since  the  secular  terms  in  ft,  w  and  o  have  argument  v 
rather  than  t  ,  they  are  not  immediately  separable  from  the  short-periodic 
terms,  whereas  for  M  the  separation  of  AM  into  Mt  and  6M  is  immediate. 


The  formulae,  in  which  arbitrary  'constants'  have  been  omitted  (see 
section  2.4)  and  in  which  compactness  is  achieved  largely  by  writing  C.  for 
cos  (jv  +  2<d)  and  for  sin  (jv  +  2u>)  ,  are:- 


6a  =  2*5  Kaq  2  |24f  +  e  [48h  cos  v  +  36f  (Cj  +  C^)] 

+  e2  [24h  cos  2v  +  18f  (2C2  +  C4)] 

+  e^  [4h  (3  cos  v  +  cos  3v)  +  3f  (C_j  +  3Cj  +  3C^  +  C,.)]  j  > 

.  (55) 


6e 

6i 

Aft 


**8K  |48h  cos  v  +  4f  (3Cj  +  70^)  +  6e  (4h  cos  2v  +  lOf  C 2  +  3f 

+  e2  [4h  (3  cos  v  +  cos  3v)  +  f  (3C_j  +  33Cj  +  17  +  3C^)]| . 


I^K  sin  2i  ^3C2  +  e  [3Cj  +  Oj]}  , 

£k  cos  i  j-6v  +  3S2  -  e  [6  sin  v  -  3S  j 


(57) 


(58) 
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Aw  ■  4*g K  e  1  -|48h  sin  v  -  4f  (3Sj  -  7Sj)  +  6e  [4  (4  -  5f)  v 

+  4h  sin  2v  -  2  (2  -  5f)  S2  +  3f  Sj 

r\ 

+  e  [6  (14  -  17f)  sin  v  +  4h  sin  3v  -  3f  S_j 

-  3  (8  -  15f)  Sj  -  (8  -  1 9f )  S3  +  3f  ,  (59) 

An  =  -  ^gK  e  1  q  |48h  sin  v  -  4f  (3Sj  -  7S3)  +  6e  [-8h  v 

+  4h  sin  2v  -  6f  S2  +  3f  f-60h  sin  v 

+  4h  sin  3v  -  f  (3S_,  +  51 Sj  +  13S3  -  3S5)]|  (60) 

and 

AM  =  n'  t  -  ^Re"1  q  |48h  sin  v  -  4f  (3Sj  -  7S3)  +  6e  [4h  sin  2v  +  3f  Sj 

2 

-  e4  [4h  (3  sin  v  -  sin  3v) 

+  f  (3S_j  +  15Sj  +  S3  -  3S5)]}  . 

.  (61) 

The  'mean  mean  motion',  n*  ,  that  appears  in  (61)  is  an  absolute  constant 
of  the  motion  that  is  related  to  the  energy  integral  as  explained  in  Refs  2  and  4 
(see  also  section  4.1).  It  is  given  by 

,2  ,3  t  c.*>\ 

n'  a'  =  y  ,  (62) 


where  the  relation  between  a'  and  (osculating)  a  is 


(63) 


This  relation  is  exact,  not  merely  first-order  in  J2  .  To  first  order  -  of 
equation  (114)  to  be  given  in  section  4.1  -  we  have 


a  =  a'  +  5a  +  ^K.  (a^/p)  |8h  +  3e^  (4h  +  3f  cos  2w)(  ,  (64) 

with  5a  given  by  (55),  so  that  a'  corresponds  to  a  particular  choice  of  mean 
element  I  . 

Two  other  comments  are  worth  making  about  equations  (55)  to  (61).  First, 
the  quantities  on  the  right-hand  sides  should  all  be  regarded  as  mean,  so  that 
strictly  K,  e,  v  should  be  written  throughout;  being  only  first-order  equations, 
however,  they  are  actually  correct  as  they  stand.  Secondly,  the  presence  of  the 
q-factor  in  (60)  and  (61)  makes  a  simple  complete  expression  for  SU  impossible 


( of  the  corresponding  remark  in  section  2.3);  complete  expressions  for  6£  and 
6n  are  simple  enough,  however,  but  will  not  be  given  in  this  paper. 

We  now  write  down  the  perturbations,  to  0(e)  only,  that  we  shall  require 

• 

in  the  sequel.  For  each  element  X,  we  express  C  and  6c  in  the  form  K£  and 
^  •  •  * 

KCj  respectively,  except  that  we  cannot  define  Mj  or  Uj  in  this  way,  as  M 

and  U  are  not  constant  for  an  unperturbed  orbit.  It  is  convenient  to  drop  e, 

w,  a  and  M  in  favour  of  ?,  n,  U  and  n  ,  and  to  introduce  k-constants  at 

appropriate  points.  Then  the  formulae  required  are:- 


al 

= 

$a  ^2(f  cos  2u  +  k^)  +  e  [4K  cos  v  +  3?  cos  (u  +  S3)  +  3f  cos  (2u  + 

•  •  •  < 

...  (65) 

*1 

= 

iV  sin  2i^3(cos  2u  +  k^)  +  e  (3  cos  (u  +  S3)  +  cos  (2u  +  v)]|  , 

(66) 

= 

-  n  cos  I  , 

(67) 

= 

£  cos  i  ^3  (sin  2u  +  k^)  -  e  [l8  sin  v  -  3  sin  (u  +  S3)  -  sin  (2u  + 

v)]3  , 

•  •  •  • 

..  (68) 

A 

5I 

= 

-  Jne  (4  -  5f)  sin  S3  , 

(69) 

51 

= 

2*5  |2  [3  (4E  +  f)  cos  u  +  7f  cos  3u]  +  3e  [6  (1  -  f)  cos  (u  +  v) 

-  2  (1  -  5f)  cos  (2u  +  S3)  +  3f  cos  (3u  +  v)  +  k^}|  , 

(70) 

• 

= 

Jne  (4  -  5f)  cos  S3  , 

(71) 

nl 

= 

2V  |2  [3  (4h  -  f)  sin  u  +  7f  sin  3u3  +  3e  [2  (1  -  3f)  sin  (u  +  v) 

-  2  (1  -  5f)  sin  (2u  +  S3)  +  3f  sin  (3u  +  v)  +  k^3| 

9 

•  •  •  • 

••  (72) 

u  < 

= 

• 

n'  +  S3  ,  ie  n’  +  ^Kn  (4  -  5f)  , 

(73) 

u> 

= 

—  2^  ^ 6  (2  -  5f)  sin  2u  -  e  [6  (26  -  33f)  sin  v 

and 

-  3  (4  -  9f)  sin  (u  +  S3)  -  (4  -  1 7 f )  sin  (2u  +  v)}^ 

(74) 

n 

1 

= 

-  |n  ^2  (f  cos  2u  +  kn)  +  e  (4h  cos  v  +  3f  cos  (u  +  S3)  +  3f  cos  (2u  + 

V)]}  . 

•  •  •  • 

..  (75) 

NB.  These  formulae  are  only  valid  to  0(e)  ;  also,  kfl  is  redefined  in  section  8. 
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The  coefficient  of  sin  v  in  (68)  is  three  times  as  large  as  in  (58),  it 
will  be  noted.  The  difference  is  due  to  the  e-term  in  (49),  which  comes  in  when 
the  v-term  in  (58)  is  replaced  by  an  M-term  that  (because  it  is  linear  with 
time)  can  then  be  removed  from  dfl  and  represented  by  0  ,  given  by  (67) . 

The  quantity  ii  that  appears  in  (67),  (69),  (71)  and  (75)  is  a  mean  element 
that,  like  all  the  mean  elements,  is  arbitrary  to  the  extent  of  a  constant,  viz 
the  constant  kn  that  appears  in  (75).  Since 

M  =  MQ  +  n'  t  , 

from  (61),  ie  n*  =  M  ,  it  may  seem  natural  just  to  take  n  »  n'  ,  but  alterna¬ 
tive  assumptions  may  be  preferable  and  will  be  considered  at  the  end  of  this 
section. 

Formulae  (69)  to  (72)  are  equivalent  to  formulae  for  the  elements  e  and 
u)  ,  viz 


Si  = 

0  , 

(76) 

el  = 

h 

|2  [l2h  cos  v  +  3f  cos 

(u  +  3;  +  7f 

cos  (2u  +  v)] 

+  3e  £<4h 

cos  2v  +  JOf 

cos  2u  +  3f  cos  2(u  +  v)  +  k£] 

}■ 

(77) 

51  = 

in 

(4  -  5f) 

(78) 

and 

“l  = 

^2e  [l2h  sin  v  -  3f  sin  (u  +  ui)  + 

7f  sin  (2u  +  v)] 

+  3  [4h  sin  2v  -  2  (2 

-  5f)  sin  2ii 

+  3f  sin  2(u  +  v)  +  » 

(79) 

where 

k 

k„  cos  ui  + 

k  sin  u 

(80) 

e 

K 

n 

and 

k  = 

-  k  sin  a>  ■ 

h  k  cos  oi 

(81) 

u 

S 

n 

It  is  because  of  the  e  '-factor  in  (79)  that  the  non-singular  and  r)|  should 

be  used  instead  of  ej  and  Wj  ,  but  tu  is  free  of  this  factor  -  in  other  words, 

u  only  has  an  O(Ke)  effect,  apart  from  its  action  in  combination  with  n'  in 

(73),  as  (69)  and  (71)  confirm.  In  fact  the  most  satisfactory  way  to  generate 

(osculating)  e  and  u  from  (mean  at  epoch)  e  and  SL  is  by  combining 
•  •  u 

formulae  for  Ij,  5  ,  and  n ^  ,  such  that 
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5  =  eQ  cos  (5q  +  cot)  +  K5j  (82) 

and 

n  =  eQ  sin  (o>0  +  Sit)  +  Kn j  »  (83) 

from  which  e  and  u>  follow  by  inversion  of  (1)  and  (2). 

We  can  now  generate  perturbation  formulae  for  quantities  more  closely 
related  to  actual  satellite  position.  It  is  convenient  to  express  these  formulae 
by  extending  the  suffix-1  notation  so  far  only  used  for  orbital  elements.  We 
start  with  r  and,  from  (40)  and  (37),  can  write 

2  2  3 

r  =  a  (1  -  e  cos  M  +  e  sin  M)  +  0(e  )  , 

ry 

r  =  a  [l  -  (£  cos  U  +  n  sin  U)  +  (£  sin  U  -  n  cos  U)  ]  +  0(e  )  .  (84) 


ze 


From  (84)  it  then  follows,  as  ir  section  6.1  which  we  are  anticipating,  that  to 
0(e)  , 


Tj  =  (r/a)  aj  -  a  (£j  cos  u  +  n j  sin  u)  +  aeUj  sin  v  , 

a  formula  equivalent  to  one  given  in  Ref  27,  viz 

fir  =  fia  (1  -  e  cos  v)  -  a  fie  cos  v  +  ae  6M  sin  v 
Now  (85)  reduces  to 

Tj  *  &a  jf  cos  2u  -  6 (h  -  kfl)  +  |e  [l2(3h  -  2k&)  cos  v 

+  9f  cos  (u  +  w)  -  3k^  cos  u  -  3k^  sin  u]|  , 

and  this  can  be  further  reduced  to  the  very  simple  formula 

r,  =  ilf  cos  2u  , 

by  making  the  following  choices  for  ka>  k^  and  k^: 

k  =  h  , 
a 


and 


kj.  =  (4h  +  3f)  cos  u 


k  =  (4h  -  3f)  sin  to 

n 


(85) 


(86) 


(87) 

(88) 

(8°) 

(90) 

(91) 


Thus  (89),  (90)  and  (9 1 )  give  'preferred  values'  for  these  three  constants  and 
they  will  be  referred  to  as  such  throughout  this  paper. 
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A  few  comments  are  worth  adding  before  we  pass  from  the  perturbation  in  r 
to  that  in  u  .  First,  the  value  of  a  that  corresponds  to  the  choice  of  k 

& 

given  by  (89)  is  often  known  as  the  'Kozai  semi -major  axis',  from  its  use  in 
Ref  6.  Second,  this  choice  of  k  is  certainly  not  the  only  useful  one: 

cl 

comparison  of  (65)  with  (64)  indicates  that  a  can  be  identified,  to  0(e)  , 
with  the  exact  constant  a'  ,  by  taking  k  =  yh  ;  some  authors  also  find  it 

Si 


convenient  to 

take 

k  =  0 
a 

.  Third,  values  of  kg  and  k^  corresponding 

to 

the  choice  of 

kc 

and  k 

n 

by 

(90)  and  (91)  are  given  at  once  by  (80)  and 

(81), 

viz 

k 

e 

=  4h  +  3f  cos  2 So 

(92) 

and 

k 

u 

=  -  3?  sin  2S 

(93) 

Fourth,  the  constants  k_,  k  ,  k  and  k  are  not  really  'constant'  at  all, 

£  n  e  a) 

since  5  has  a  secular  variation,  but  k^  and  k^  are  multiplied  by  e  in 

(70)  and  (72),  so  their  variation  involves  0(Ke)  terms  which  we  are  neglecting 

(see  also  section  9).  Finally  (and  this  ties  up  the  last  comment  with  an  earlier 

remark),  it  is  sometimes  convenient  (eg  in  Ref  27)  to  suppress  the  secular 

•  JL  i 

variation  5  ,  even  in  (82>  and  (83),  using  for  M  a  quantity  that  is  really  U  ; 

this  involves  the  use  of  an  incorrect  value  of  f?  ,  and  hence  of  r  (from  the 

standard  algorithm  of  section  2.5),  but  this  can  be  allowed  for  by  addition  of 

an  extra  term,  viz  - jient  (4  -  5f)  sin  v  ,  to  (87). 

We  now  turn  to  u  ,  and  from  (49)  can  write 

u  =  U  +  2e  sin  M  +  lj  e2  sin  2M  +  0(e3)  , 

ie 

u  =  U  +  2  ((  sin  U  -  n  cos  U)  +  |J  [(£2  -  n2)  sin  2U  -  2£n  cos  2u]  +  0(e3)  . 

.  (94) 


This  gives,  correct  to  0(e) 
Uj  ■  Uj  +  2  (£,  sin  u  -  n j 


,  again  anticipating  section  6.1, 

u)  +  ie  ^4U j  cos  v  +  Cj  [sin  (u  +  v)  +4  sin  So] 

-  rij  [cos  (u  +  v)  +4  cos  5)]| 


cos 


(95) 


which  reduces  to 
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Uj  *  -  i*j  j(6  -  7f)  sin  2u  -  I  [6(8  -  9f)  sin  v  -  (6  -  f)  sin  (u  +  tH) 

-  2(1  -  f)  sin  (2u  +  v)  +  3k^  sin  u 

~  3^  cos  u]|  .  (96) 

On  substituting  the  preferred  values  of  k^  and  k^  ,  given  by  (90)  and 
(91),  we  obtain 

Uj  =  -  iV  |(6  -  7f)  sin  2u  -  2e  [6(5  -  6f)  sin  v  -  (3  -  5f)  sin  (u  +  w) 

-  (1  -  f)  sin  (2u  +  v)]}  .  (97) 

If  iu  is  suppressed  in  (82)  and  (83),  an  extra  term  is  needed,  viz 
-ent  (4  -  5f)  cos  v  . 

We  are  now  equipped  to  derive  the  perturbations  in  the  cylindrical  polar 

coordinates  (r',  u' ,  c)  introduced  in  section  2.5.  From  (43),  (66)  and  (68), 

we  have  that,  to  0(e)  , 

Cj  =  -  ir  sin  2i  |(1  -  k^)  sin  u  +  kfi  cos  u  -  l^e  [2  sin  (u  +  v)  -  3  sin  . 


If  we  make  the  preferred  choice 


k.  =  1 

l 


and  also  kfi  =  0  (it  is  arguable  that  k^  should  not  have  been  introduced 

anyway,  since  it  should  naturally  operate  as  a  ’coefficient  of  sin  0u'  whereas 

k  and  k.  operate  as  coefficients  of  cos  Ou  ),  (98)  reduces  to 
a  i 


c.  =  -$■?£  sin  2i  [2  sin  (u  +  v)  -  3  sin  w] 


(100) 


Now  we  know  from  (46)  that  we  take 


rj  =  ri  , 


(101) 


so  it  remains  to  deal  with  u’  .  We  define  uj  on  the  basis  that  u'  and  u 
are  identical;  then  (97)  (44)  and  (68)  give,  with  the  preferred  constants, 


uj  =  rS  | f  sin  2u  +  4e  [6h  sin  v  +  f  sin  (u  +  w)]|  , 

and  this  could  also  have  been  obtained  from  the  formula,  correct  to  0(e)  , 


002) 


1  uu 
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6U  +  Se  (2  sin  v  +  Je  sin  2v)  +  e  6M  [2  cos  v  +  e  (2  +  $  cos  2v)] 


+  60  cos  i  ,  (103) 

which  is  equivalent  to  the  first  formula  on  page  36  of  Ref  27  . 

From  uj  may  be  derived  the  quantity  associated  with  the  'L' 

coordinate  referred  to  at  the  end  of  section  2.5,  though  we  have  seen  that  the 
cartesian  'RLC  system’  of  coordinates  is  less  satisfactory  than  the  cylindrical 
'RUC'  system  given  by  r' ,  u'  and  c  .  Thus  (47)  gives 


=  j*t  r  |f  sir.  2u  +  4e  [6h  sin  v  +  f  sin  (u  +  ii>)jV  , 


(104) 


it  being  irrelevant,  to  first  order,  whether  (102)  is  multiplied  by  r  or  r  ; 
on  expressing  in  terms  of  a  (rather  than  r  or  r  ),  however,  we  get 

=  5*3  a  |2f  sin  2u  +  e  f48h  sin  v  +  7f  sin  (u  +  S3)  -  f  sin  (2u  +  v)J j 

(In  (100),  by  contrast,  r  can  be  replaced  by  either  r  or  5.) 

Preferred  values  have  now  been  given  for  all  the  constants  except  , 

introduced  in  (75).  The  value  of  depends  on  whether  or  not  it  is  regarded 

as  essential  to  preserve  the  familiar  Kepler-third-law  relation,  ie  to  demand 
-2  -3 

that  n  a  =  y  (of  (4)  and  (62)).  If  we  do  want  to  preserve  this  relation,  then 
it  follows  at  once  from  (65)  and  (75)  that  we  must  set  k^  equal  to  k^  .  If 
we  want  to  assign  ka  and  k^  independently,  on  the  other  hand,  then  the  Kepler 
relation  must  be  replaced  by  the  more  general  formula 


-2  -3 
n  a 


„  [l  -  3r  <ka  -  v] 


(105) 


It  has  been  remarked  that  "it  may  seem  natural  just  to  take  n  =  n'  " 
(since  n'  *  M  ),  and  we  have  also  seen  that 


k.  -45 


from  which  it  follows  that 


n  *  n'  =*  kR  =  $-h 

But  h  ,  not  i’h  ,  is  the  preferred  value  of  k  ,  corresponding  to  the  Kozai 

semi-major  axis,  a^  say.  This  would  lead  to  a  'mean  mean  motion'  of  no 

_____  .  •  2 
*  This  follows  on  replacing  6L  in  Ref  27  by  6U  +  60  cos  i  +  |e  6M  ;  the  last 

term  is  0(e)  and  should  have  been  included  in  the  expression  for  6L  on 

page  2  of  Ref  27. 
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practical  use,  so  it  is  common  practice*  to  combine  k  -  h  with  k  «  .  in 

a  n 

which  case 


-2  -3  ,23  -r. 

na  =  n'  a  -  n  (1  -  Kh)  . 


The  general  relation  connecting  n,  k^  and  n'  is 


n  =  n'  +  Kn  (1J  kR  -  h) 


006) 


(107) 


and  we  will  use  this  to  derive  a  value  for  k^  that  for  much  of  this  paper  is  not 
merely  'preferred'  but  actually  essential.  The  associated  value  of  n  is  given 
by 


n  =  U  (=  n'+uby  (73))  , 


(108) 


the  significance  of  this  value  being  that  it  permits  the  second-order  development 
of  section  4  to  proceed  on  the  basis  that 

2n  J  (cos  2U,  sin  2U)  dt 

can  be  replaced  by  (sin  2U,  -cos  2U) ,  apart  from  the  matter  of  integration  con- 
stants.  Since  w  is  given  by  (78),  (107)  now  yields 


kn  =  (3  -  4f )  , 


from  which  (105)  gives 


n2  a^  *  v  [l  +  JK  (6  -  7f)]  . 


(109) 


(110) 


If,  on  the  other  hand,  is  not  limited  to  the  value  given  by  (109),  then  (73), 

(78)  and  (107)  lead  to  the  general  relation 


U  -  n  =  jKn  (6  -  8f  -  3kR) 


(111) 


We  may  regard  the  right-hand  side,  here,  as  defining  KUj  ,  since  Uj  was  not 
defined  when  the  other  were  introduced. 


*  76  —  -  - - - - 

Kozai  effectively  adopts  a  set  of  k's  as  follows:  ka  and  ke  are  as  given 

by  (263)  and  (264)  of  section  8,  such  that  ka  agrees  with  our  preferred  value 
(See  section  8)  to  0(e)  but  ke  has  an  'h  component’  that,  to  0(e)  ,  is  double 
our  value  given  by  (92);  k£  is  0  (of  1 ,  here);  kn  is  0  ,  as  here;  k^  is 
exactly  as  given  by  (93)  here;  finally  kM  (not  defined  till  section  8)  is 
identical  with  k„,  . 
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4  SECOND-ORDER  PERTURBATIONS  IN  OSCULATING  ELEMENTS 

For  any  element,  t,  ,  the  secular  rate  of  change  may  be  expressed  in  the 


oo 

*  ■  £ 


(112) 


and  the  short-periodic  perturbation  in  the  form 


uu 

■  £  ^  <5  ' 


(113) 


generalizing  the  notation  of  section  3.  We  are  now  concerned  with  formulae  for 
^2  and  t, ^  ,  but  will  be  neglecting  0(e)  contributions. 

The  basic  idea  is  to  'bootstrap*  on  the  first-order  solution,  substituting 
it  on  the  right-hand  sides  of  the  planetary  equations  and  re-integrating.  The 
planetary  equations  must  be  taken  to  0(e)  in  this  process,  since 

K  e  =  K  (e  +  6e)  , 


where 


<5e  =  K  ej  +  0(K  )  , 


and  hence  Ke  terms  lead  to  K  terms  without  the  factor  e  .  A  special 
procedure,  free  of  further  integration,  is  possible  for  £  =  a  ,  and  this  is 
developed  first.  It  may  be  compared  with  a  re-derivation  of  a2  by  the  general 
procedure. 

4 . 1  Perturbation  in  a  (special  method) 

From  (63)  and  (11)  it  follows  that 


a  =  a ' 


+  K  a'  q  ^  (p/r)^  (f  cos  2u  +  ^h)  , 


(114) 


this  being  an  exact  result.  Recovery  of  the  first-order  a^  ,  given  by  (65)  to 
0(e)  ,  is  immediate  on  expanding  p/r  by  (12),  on  the  basis  that 


=  r.  (1  +  -  fh)  . 


(115) 


If  we  are  to  be  able  to  write 


=  a'  +  K  5  q  ^  (p/r)^  (f  cos  2u  +  |h)  +  £^2  , 
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there  will  be  five  sources  of  terms  in  :  (i)  the  'a  variation',  lost  on 

replacing  K  by  K  ;  (ii)  the  generalization  from  a'  to  I  ,  these  being 
connected  by  (115);  (iii)  the  perturbation  in  (p/r)  ;  (iv)  the  variation  in  i  , 
which  affects  f  and  h  ;  and  (v)  the  perturbation  in  u  .  The  contributions 
from  these  sources  are  as  follows,  neglecting  terms  that  are  0(e)  in  a„ 


(i)  -  2a  (f  cos  2u  +  $h)  (f  cos  2u  +  k  )  ; 

a 

(ii)  a  (f  cos  2u  +  $h)  (k  -  3h)  ; 

d 

(iii)  3a  (f  cos  2u  +  ih)  (h  +|f  cos  2u)  ; 

(iv)  If  (1  -  f  )  (cos  2u  -  1 ) (cos  2u  +  k^)  ; 

(v)  £lf  (6  -  7f )  sin2  2u  . 


Summation  of  the  contributions  gives 

a^  =  ia  |f2  cos  4u  +  f  [5  -  9f  -  3k&  +  3k^(l  -  f)]  cos  2u 

+  4  [14  -  33f  +  24f 2  -  6kah  -  9k.f(l  -  f )]  }  .  (116) 

The  constant  term  here  is  of  no  great  general  importance,  and  in  the  next  section 
an  arbitrary  constant  k0  is  introduced,  by  equation  (121),  in  analogy  with  k  . 
Comparison  of  (116)  and  (121)  enables  us  to  obtain  the  expression  for  k£a 
appropriate  to  an  interpretation  of  I  as  a'  ,  ie  to  obtain  the  second-order 
relation  of  a  to  a'  .  We  set  k&  to  §h  ,  therefore,  and  obtain 

k2a  =  *  C(l°  “  21*  +  I5f2)  "  9k.f(l  ~  £>]  .  (117) 

An  application  for  this  result  is  found  in  section  9.7. 

4.2  Perturbation  in  a  (general  method) 

The  starting  point  is  the  exact  equation  (23).  From  this  it  follows,  using 
(12),  (51)  and  (52),  that  to  0(e) 

a  =  -  jKna  ^4f  sin  2U  +  e  [4h  sin  v  -  f  sin(u  +  w)  +  21 f  sin  (2u  +  v)]| 

.  (118) 

The  object  of  replacing  sin  2u  by  sin  2U  in  the  main  term  is  to  permit  first- 
order  integration  with  respect  to  time  (as  opposed  to  changing  the  integration 
variable  to  true  anomaly,  by  invoking  (54)).  The  first-order  solution  then  follows 
at  once  in  the  form 
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6a 


=  $Ka  ^ 2f  cos  2U  +  e  [4h  cos  v  -  f  cos  (u  +  C)  +  7f  cos  (2u  +  v)]|  , 


(119) 


which  is  equivalent  to  (65),  with  kfl  ■  0  ,  if  allowance  is  made,  by  (53),  for 
the  difference  between  cos  2U  and  cos  2u  . 

It  is  now  possible  to  expand  A  ,  as  given  by  (118),  in  terms  of  the  first- 

order  reference  solutions  for  a,  e,  i,  u>,  U  and  n  .  Thus  there  are  six 

contributions  to  the  total  A  which  are  second-order  and  they  may  be  de.ioted  by 

A  a,  A  a  etc,  where 
a  e 

*  •  3a  Ar 

V  '  ■ 

2 

Remembering  that  K  is  a  function  of  a  ,  and  that  we  are  neglecting  0(K  e) 
perturbations,  we  obtain 


and 


A  a 
a 


A  A 
e 


A. a 

l 


A  a 
01 


A  a 
U 


A  a 
n 


-  Ka  (f  cos  2u  +  k  )  , 

a 


2 

=  -  na  [4h  sin  v  -  f  sin  (u  +  o>  )  +  21  f  sin  (2u  +  v)]  x 
x  (j2h  cos  v  +  3f  cos  (u  +  to)  +  7f  cos  (2u  +  v)]  , 

=  Ka  (1  -  f ) (cos  2u  +  k^)  , 

2 

=  jVk  na  [4h  cos  v  +  f  cos  (u  +  oi)  +  21  f  cos  (2u  +  v)J  x 

x  [l  2h  sin  v  -  3f  sin  (u  +  ui)  +  7f  sin  (2u  +  v)J  , 

=  -  iKa  (2  -  5f )  cos  2u 


-  IjKA  (f  cos  2u  +  k  ) 


It  is  immaterial  whether  the  quantities  on  the  right-hand  sides  are  regarded  as 
osculating  or  mean,  so  'bars'  are  omitted  for  convenience.  The  value  of  kR  must 
be  taken  from  (109),  however. 

It  now  follows,  after  some  tedious  reduction,  that  the  total  second-order 
contribution  to  a  is  given  by 

Aa  =  -$K2naf  ^2f  sin  4u  +  [5  -  9f  -  3kfl  +  3k£(l  -  f)]  sin  2u|  .  (120) 


o 

o 


The  formula  for  second-order  6a  follows  immediately,  and  we  may  write 


a2  -  4a  {f2  cos  4u  +  f  [5  -  9f  -  3k&  +  3k.(l  -  f>]  cos  2u  +  3k2a|  , 

.  (121) 

where  it  has  been  convenient  to  introduce  a  second-order  constant,  k2a  ,  akin 
to  k  .  Clearly  (121)  is  in  agreement  with  (116). 

SL 

Taking  the  preferred  values  of  h  and  1  for  kfl  and  k^  respectively, 

we  get 


a2  =  4a  (f2  cos  4u  +  5fh  cos  2u  +  3k2a)  ,  (122) 

where  the  preferred  value  of  k2a  ,  as  we  shall  see  in  section  6.1,  is  given  by 

k,  -  h2  -  7*5  f 2  ,  (123) 

za 

ie  by 

k2a  =  1  -  3f  +  244f2  . 

4.3  Perturbation  in  i 

The  starting  point  is  the  exact  equation  (25),  which  to  0(e)  gives 

i  =  -  $Kn  sin  2i  j 2  sin  2U  -  e  [sin  (u  +  to)  -  7  sin  (2u  +  v)]|  .  (124) 

The  first-order  solution  follows  at  once,  in  the  form 

6i  =  T^K  sin  2i  |3  cos  2U  -  e  [3  cos  (u  +  w)  -  7  cos  (2u  +  v)]|  ,  (125) 

equivalent  to  (66)  (with  k^  =  0)  . 

* 

We  expand  i  ,  as  given  by  (124),  in  terms  of  the  first-order  reference 
solutions  for  a,  e  eta,  just  as  in  section  4.2,  obtaining  six  second-order 
contributions  to  Ai  ,  viz 
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A  i  =  -  2Ki  (f  cos  2u  +  k  )  , 

ci  d 

Agi  =  n  sin  2i  [sin  (u  +  w)  -  7  sin  (2u  +  v)}  x 

x  [l2h  cos  v  +  3f  cos  (u  +  u>)  +  7f  cos  (2u  +  v)}  , 

•  • 

A^i  =  ^Ki  cos  2i  (cos  2u  +  k^)  , 

*  2  r* 

A  i  =  n  sin  2i  [cos  (u  +  id)  +7  cos  (2u  +  v)l  x 

CO  u  J 

x  [l2h  sin  v  -  3f  sin  (u  +  u>)  +  7f  sin  (2u  +  v)J  , 

Ayi  =  -  ^Ki  (2  -  5f)  cos  2u 

A  i  -  -  l^Ki  (f  cos  2u  +  k  ) 
n  n 

Tedious  reduction,  as  for  Aa  ,  leads  to  the  total  Ai  given  by 

Ai  «  25K^n  sin  2i  j(3  +  5f)  sin  4u  +  6  [f  +  4ka  -  k^(l  -  2f)J  sin  2u|  , 

.  (126) 

from  which  the  formula  foi  second-order  6i  follows  immediately;  thus 

i2  =  ~  5*2  sin  2i  |(3  +  5f)  cos  4u  +  12  [f  +  4ka  -  k^,(l  -  2f)]  cos  2u  +  k^j  , 

.  (127) 


and 


with  a  convenient  second-order  constant  introduced. 

With  the  preferred  values  of  k  and  k.  ,  we  get 

a  1 

i2  =  -  9*2  sin  2i  [(3  +  5f)  cos  4u  +  36(1  -  f)  cos  2u  +  k^]  ,  (128) 

where  the  preferred  value  of  k^  is  given  (as  we  shall  see  in  section  6.3)  by 


k,.  =  33  -  49f 

2i 


(129) 


4.4  Perturbation  in  Q 


The  exact  starting  equation  is  (26),  which  to  0(e)  gives 

U  -  -  jKn  cos  i  ^2  (I  -  cos  2U)  +  e  [6  cos  v  +  cos  (u  +  u>)  -  7  cos  (2u  +  v)j|  . 

.  (130) 
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The  first-order  solution  follows  at  once,  and  is  composed  of  a  secular  term, 
with  given  by  (67),  and  a  periodic  term  in  the  form 


5fi 


cos  i  (3  sin  2U  -  e 


[18 


sin  v 


+  3  sin  (u  +  to)  -  7  sin  (2u 


*>]} 


031) 


equivalent  to  (68)  (with  =  0) . 

On  expanding  $  relative  to  the  first-order  reference  solutions  for  a, 
e  eta ,  we  obtain  the  following  six  second-order  contributions  to  A&  : 

A  H  =  -  2Kft  (f  cos  2u  +  k  )  , 

a  a 

2 

n  cos  i  [6  cos  v  +  cos  (i>  +  to)  -  7  cos  (2u  +  v)}  x 

x  [l2h  cos  v  +  3f  cos  (u  +  to)  +  7f  cos  (2u  +  v)]  , 

A^ft  =  -  jKftf  (cos  2u  +  k^)  , 

#  2 

A^ft  =  -  n  cos  i  [6  sin  v  -  sin  (u  +  to)  -  7  sin  (2u  +  v)]  * 

x  [l2h  sin  v  -  3f  sin  (u  +  to)  +  7f  sin  (2u  +  v)]  , 

A^ft  ■  -  Kft  (2  -  5f)  cos2u 

and 

A  ft  =  -  1 |Kft  (f  cos  2u  +  k  ) 

n  2  n 

The  total  Aft  is  given  by 

Aft  »  -  iK^n  cosi  j(3  +  f)  cos  4u  -  3(2f  -  4k^  -  k^f)  cos  2u 

+  (15  -  1 9f  -  12k  -  3k. f  -  9k  )1  ; 

a  l  n  / 

.  0  32) 

thus  there  is  a  secular  component  of  the  second-order  Aft  ,  and  kR  has  been 
left  in  (132)  because  the  value  required  depends  on  the  chosen  n  ,  which  appears 
explicitly  in  ft  .  We  get 

• 

ft.  =  -in  cos  I  (15  -  19f  -  12k  -  3k. f  -  9k  )  (133) 

2  a  i  n 

and 

ft.  ■  -  ’h  cos  i  [(3  +  f)  sin  4u  -  6(2f  -  4k  -  k.f)  sin  2u]  ,  (134) 

i  a  i 

no  k.  constant  being  required  (it  would  be  the  coefficient  of  sin  Ou  !) 

2ft 


With  the  preferred  values  of  kfl  and  k^  , 


we  get 


(135) 


^2  “  -  £n  cos  i  (3  -  4f  -  9k^) 

and 

q  =  -  I1?  cos  I  [(3  +  f)  sin  4u  +  6(4  -  7f)  sin  2u]  .  (136) 

With  our  usual  value  of  k^  ,  we  get  (of  (271)  in  section  9.5) 

^2  "  ■§"  n  cos  i  (3  -  4f)  .  (137) 

4.5  Perturbation  in  g 

The  exact  starting  equation  is  (31),  which  to  0(e)  gives 

l  m  -  jKn  |(4h  +  f)  sin  U  +  7f  sin  3U 

+  2e  [(4  -  5f)  sin  u>  +  (7  -  8f)  sin  (u  +  v) 

-  (1  +  2f )  sin  (2u  +  ui)  +  17f  sin  (3u  +  v)]J  . 

.  (138) 

The  first-order  solution  follows  at  once,  in  a  form  consisting  of  a  secular  term, 
with  t]  given  by  (69) ,  and  a  periodic  term 

6g  -  ^6  (4h  +  f)  cos  U  +  14f  cos  3U  +  3e  x 

x  [2(7  -  8f )  cos  (G  +  v)  -  2(1  +  2f)  cos  (20  +  u)  +  17f  cos  (3u  +  v)]} 

.  (139) 

equivalent  to  (70)  with  k^  taken  as*  2(4h  +  f)  cos  u  . 

On  expanding  t  relative  to  the  usual  first-order  reference  solutions, 
we  obtain: 


This  k  and  the  corresponding  k^  (see  section  4.6)  are  equivalent  to  the 
(effective)  kg  and  ky  of  Ref  25  that  we  encounter  in  section  5.2. 
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V 


A  £ 

CO 


and 


V 


-  2K£  (f  cos  2u  +  k  )  , 

8L 


-  n  [(4  -  5f)  sin  o>  +  (7  -  8f)  sin  (u  +  v)  -  (1  +  2f)  sin  (2u  +  oj) 

+  17f  9in  (3u  +  v)]  [l2h  cos  v  +  3f  cos  (u  +  to)  +  7f  cos  (2u  +  v)]  , 

$K2nf  (1  -  f ) (5  sin  u  -  7  sin  3u) (cos  2u  +  k^)  , 

-  Ak  n  [(4  -  5f)  cos  to  -  (7  -  8f)  cos  (u  +  v)  -  (1  +  2f)  cos  (2u  +  to) 

-  17f  cos  (3u  +  v)][l2h  sin  v  -  3f  sin  (u  +  to)  +  7f  sin  (2u  +  v)]  , 

I^K2n  (2  -  5f)  [(4h  +  f)  cos  u  +  2 1  £  cos  3u]  sin  2u 


-  1$K£  (f  cos  2u  +  kR) 


The  total  A£  is  given  by 

A£  -  ^-6K2n  ^5f  (14  -  17f)  sin  5u  +  [72  -  346f  +  449f2  +  336fk  -  168k.f(l  -  f)]  * 

l  a  i 

X  sin  3u  -  2[l08  -  200f  +  77f2  -  24kfl(4  -  5f)  -  60^(1  -  f)]  sin  u|  , 

.  O40) 

whence 

C2  -  ~  2+8  {3f  (14  -  1 7f )  cos  5G  +  [72  -  346f  +  449f2  +  336fka 

-  168k. f(l  -  f)]  cos  3u  -  6 [l08  -  200f  +  77f2  -  24k  (4  -  5f) 
x  a 

-  60k^f(l  -  f )]  cos  u  |  , 

.  041) 

no  constant  being  introduced. 

With  the  preferred  values  of  k  and  k.  ,  we  get 

a  x 


-  rta  { 


3f  (14  -  1 7f)  cos  5u  +  (72  -  178f  +  1 1 3f  )  cos  3u 

-  6(12  +  4f  -  43f2)  cos 


=  } 


042) 


mi 
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The  first-order  solution  follows  at  once,  in  a  form  consisting  of  a  secular  term, 
with  rij  given  by  (71),  and  a  periodic  term 

6n  =  |6(4h  -  f)  sin  U  +  14f  sin  3U  +  3e  x 

x  [10(1  -  2f )  sin  (u  +  v)  -  2(1  +  2f)  sin  (2u  +  5)  +  17f  sin  (3u  +  v)]|  , 

.  (144) 

equivalent  to  (72)  with  k^  taken  as  2(4h  -  f)  sin  5  . 

On  expanding  fi  relative  to  the  usual  first-order  reference  solutions,  we 
obtain: 

An53-  2Kti  (f  cos  2u  +  k  )  , 

cl  3 

Agn  =  ^?K^n  [(4  -  5f)  cos  u  +  5(1  -  2f)  cos  (u  +  v)  -  (1  +  2f)  cos  (2u  +  io) 

+  17f  cos  (3u  +  v)]  [l2h  cos  v  +  3f  cos  (u  +  io)  +  7f  cos  (2u  +  v)]  , 

o 

A^n  ■  -  1|K  nf  (1  -  f)  (cos  u  -  cos  3u)(cos  2u  +  k^)  , 

Awn  =  -  5*?K2n  [(4  -  5f)  sin  ui  -  5(1  -  2f)  sin  (u  +  v)  -  (1  +  2f)  sin  (2u  +  u) 

-  1 7 f  sin  (3u  +  v)][l2h  sin  v  -  3f  sin  (u  +  w)  +  7f  sin  (2u  +  v)]  , 

Ayfi  =  (2  -  5f )  [(4h  -  f)  sin  u  +  21f  sin  3u]  sin  2u 

and 

A  n  ■  -  1  $Kn  (f  cos  2u  +  k  ) 

n  *  n 

The  total  An  is  given  by 

An  -  -  ?VK2n  {5f  (14  -  1 7f)  cos  5u  +  [72  -  310f  +  395f2  +  336fka 

-  168kif(l  -  f)]  cos  3u  -  2 [84  -  224f  +  155f2 

-  24ka(4  -  7f)  -  84kif(l  -  f)]  cos  uj  ,  (145) 

whence 

n2  -  -  Th  ^3f  (14  -  17f)  sin  5u  +  [72  -  310f  +  395f2  +  336fka 

-  168k.  f(l  -  f)l  sin  3u  -  6[84  -  224f  +  155f2 
i 

-  24k  (4  -  7f )  -  84k. f(l  -  f)l  sin  u}  ,  (146) 

a  17 

no  constant  being  introduced. 

D 
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With  the  preferred  values  of  k  and  k.  ,  we  get 

n2  ■  ”  rh  {3f  (14  -  1 7f )  sin  5u  +  (72  -  142f  +  59f2)  sin  3u 

+  6(12  -  4f  +  13f2)  sin  u  |  .  (147) 

4.7  Perturbation  in  (a  +  to) 

From  the  exact  equations  (27)  and  (28)  it  follows  that,  to  0(e)  , 

a  +  to  *  i"Kn  |8  [(3  -  4f)  -  (1  -  4f)  cos  2U]  +  e  [2(38  -  5 1  f )  cos  v 

+  (4  -  17f)  cos  (u  +  to)  -  7(4  -  17f)  cos  (2u  +  v)j| 

.  (148) 

The  first-order  solution  follows  at  once;  it  consists  of  a  secular  term,  given 
by 

Oj  +  to j  “  n  (3  -  4f)  ,  (149) 

and  a  periodic  term  in  the  form 

6a  +  6oo  =  -  &K  |l2  (1  -  4f)  sin  2U  -  e  [6(38  -  51f)  sin  v 

+  3(4  -  17f)  sin  (u  +  5)  -  7(4  -  !7f)  sin  (2u  +  v)]J  .  (150) 

This  is  equivalent,  by  (53),  to  a  result  that  could  be  added  to  equations  (65)  to 
(75),  viz 

Oj  +  (0j  =  -  1*8  ^24  (1  -  4f)  sin  2u  -  e  [6(38  -  51  f)  sin  v 

-  3(4  -  15f)  sin  (u  +  to)  -  (4  -  23f)  sin  (2u  +  v)]|  .  (151) 

On  expanding  (a  +  to)  relative  to  the  first-order  reference  solutions  as 
usual,  we  obtain: 

A  (a  +  to)  =  -  2K  (cr  +  to)  (f  cos  2u  +  k  )  , 

3i  3 

o 

Ae(0  +  u)  -  ^K  n  [2(38  -  5 1  f )  cos  v  +  (4  -  17f)  cos  (u  +  to) 

-  7(4  -  1 7 f )  co 8  (2u  +  v)]  [l2h  cos  v  +  3f  cos  (u  +  to) 

+  7f  cos  (2u  +  v)]  , 

A^(a  +  to)  »  -  4K2nf  (1  —  f )  ( 1  -  cos  2u)(cos  2u  +  k^)  , 
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A  (o  +  <I>)  *  [2(38  -  51 f)  sin  v  -  (4  -  17f)  sin  (u  +  w) 

u 

-  7(4  -  17f)  sin  (2u  +  y)][l2h  sin  v  -  3f  sin  (u  +  «) 

+  7f  sin  (2u  +  v)}  , 

A^(cr  +  i )  ■>  -  $K2n  (1  -  4f)(2  -  5f)  sin22u 

and 

A r(o  +  i)  ■  -  1 JK  (o  +  u) (f  cos  2u  +  kn) 

The  total  A(o  +  u>)  is  given  by 

A(o  +  u)  -  ^K2n  | (24  -  4f  -  73f2)  cos  4u  -  4[f(64  -  5 1  f )  -  24ka(l  -  4f) 

-  48kif(l  -  f)]  cos  2u  +  [(432  -  I052f  +  637f2)  * 

-  24(3  -  4f)(4ka  +  3kn>  -  192^(1  -  f)]^  , 

. r  0  52) 

where  a  term  in  k^  has  been  retained  in  the  secular  component  for  the  same 
reason  as  with  An  .  We  get 

(o  +  3)  -  An  [(432  -  1052f  +  637f2)  -  24(3  -  4f ) (4kfl  +  3kR)  -  192k.f(l  -  f)] 

.  0  53) 

and 

(a  +  w)2  -  rh  {(24  -  4f  -  73f2)  sin  4u  -  8[f(64  -  51 f)  -  24ka(l  -  4f) 

-  48kif(l  -  f)]  sin  2G  }  . 

.  0  54) 

With  the  preferred  values  of  kfl  and  k^  ,  we  get 
(5  +  S)2  -  An  [(144  -  428f  +  253f2)  -  72kn(3  -  4f)]  (155) 

and 

(o  +  u))2  *  rfa  {(24  -  4f  -  73f2)  sin  4u  +  8(24  -  148f  +  147 f 2 )  sin  2uj  , 

.  (156) 

and  if  we  use  our  usual  value  of  k  we  get 

n 

(3  +  b2  ”  -  An  (288  -  724f  +  515f2)  . 
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4.8 


Perturbation  in  n  and  in 


To  be  able  to  analyse  the  perturbation  in  U  ,  as  v«e  shall  in  section  4.9, 

t 

we  need  a  formula  for  the  perturbation  in  J  n  dt  ,  to  be  able  to  combine  with 

0 

A(o  +  id)  ,  following  (20)  and  (3).  Thus  we  first  need  the  formula  for  n2  that 
corresponds  to  the  perturbation  in  n  . 

The  perturbation  in  n  could  be  obtained  by  the  usual  method,  involving 
integration,  but  it  is  more  direct  to  obtain  it  from  the  perturbation  in  a  , 
making  use  of  the  fundamental  Kepler  relation  (4) .  This  is  possible  even  though, 
because  of  choosing  independent  k-constants  for  n  and  a  ,  n  and  a  will  not 
in  general  satisfy  the  relation  themselves.  We  introduce  p^  and  p2  »  therefore, 
dependent  only  on  f  (so  long  as  we  work  only  to  first  order  in  e) ,  such  that 

52  I3  =  p  (1  +  K  Cj  +  K2  y2)  .  (158) 


This  is  equivalent  to  the  relation 

n2  a3  =  n2  I3  (1  +  K  Pj  +  K2  p2)  1  ,  (159) 

where  a  strong  reason  for  introducing  the  relation  in  this  form,  rather  than 
reciprocally,  is  that  the  final  expression  for  p2  -  see  (185)  -  becomes  very 
simple;  reciprocally,  the  relation  is 

n2  a3  =  n2  I3  [l  -  K  Pj  -  K2  (p2  -  p2)]  .  (160) 


The  general  first-order  result,  previously  given  by  (105),  may  now  be  written  in 
the  form 


"  3  (kn  “  V 


which  with  the  preferred  kfl  and  kn  reduces  ( of  (110))  to 


But 


and 


P,  =  i  (6  -  7f )  . 

_  -2 
a  =  a  +  K  a,  +  K  a 
1  2 

-2 

n  =  n+Kn^+K  n2» 


(161) 


(162) 
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from  which  it  follows  that 


n2  a3  =  n2  S3  O  +  K®.  +  K2^)3  (1  +  Kn.  +  K2A,)2  , 


063) 


on  defining  a^  to  be  a^/a  etc  for  convenience.  On  comparing  (160)  and  (163) 


it  follows  that 


-  (3a j  +  2nj) 


064) 


^ 2  A  A  a2  zv 

6a j  +  6ajiij  +  3n^  -  3a2  ~  2n 2 


(165) 


There  is  nothing  new  in  (164),  which  just  leads  to  (161)  again,  but  (165) 
leads  to  the  desired  formula  for  n2  •  To  see  this,  we  start  with  formulae  for 
aj  ,  nj  and  a 2  ,  which  follow  from  (65),  (75)  and  (121)  respectively,  where¬ 
upon  (165)  gives 


2n2  +  P2 


=  i  ^7f 2  cos  4u  -  4f  [2(5  -  9f)  -  12ka  +  6^(1  -  f)  -  9kJ  cos 


+  3  (5f2  +  16k2  -  24k  k  +  18k2  -  8k  )  X  . 

a  an  n  2a  J 


(166) 


(It  is  legitimate  to  leave  kn  in  the  coefficient  of  cos  2u  ,  in  spite  of  the 
remarks  of  section  3,  because  it  arises  through  expansion  of  (165)  and  not  from 
integration. ) 


But  y2  must  be  independent  of  u  ,  so  we  must  in  consequence  of  (166)  be 


able  to  write 


n0  =  I^n  Wf2  cos  4u  -  4f  [2(5  -  9f)  -  12k  +  6k.  (1  -  f)  -  9k  J  cos  2u  +  16k„ 

a  i  n  i 


(167) 


where  the  introduced  constant  k„  is  such  that 

2n 

5-  -  i  (5f2  +  16k2  -  24k  k  +  lSk*  -  8k  )  -  2k  . 

I  a  a  n  n  la  zn 


(168) 


Jith  the  preferred  values  of  k  ,  k.  and  k  ,  (167)  reduces  to 

a  i  n 


fVn  | 


7f  cos  4u  +  8f  (7  -  9f)  cos  2u  +  16k„ 


(169) 


and  ( 3  68)  to 


u2  =  I  (40  -  104f  +  73r  -  8k2a)  -  2k2n  . 


(170) 
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t 

To  obtain  the  perturbation  in  [  n  dt  ,  we  observe  that 

n 

where  nj  -  of  (75)  and  (119)  - 

n,  =  -  1  In  \f  cos  2U  +  k  + 

1  l  n 

.  (172) 

Then  writing  the  required  integral  as  a  sum  of  secular  and  periodic  components, 
we  have 

t 

[n  dt  =  (n  +  An)  t  ^  ic  /,  +  K2  /2  ,  (173) 

6 

where  (172)  and  (169)  imply  that 

An/n  =  -  It  ic  k  +  K2  k„  .  (174) 

n  2n 

For  /( ,  we  integrate  the  periodic  terms  of  (172),  obtaining 


0 

=  n  +  K  nj  +  K2  n2  ,  (171) 

is  given  by 

} e  [4h  cos  v  -  f  cos  (u  +  S)  +  7f  cos  (2u  +  v)]l  • 


Jj  =  -  i  ^3f  sin  2U  +  e  £l  2h  sin  v  -  3f  sin  (u  +  2)  +  7f  sin  (2u  +  v)]| 


(175) 


For  /2  ,  similarly,  we  integrate  (167)  to  obtain  (since  kR  must  necessarily 
be  given  its  preferred  value). 


J2  =  A  {?f2 

which  -  or  from  (169) 
well,  reduces  to 


sin  4u  +  16f  [4  -  3f  +  6kfl  -  3k_^(l  -  f)}  sin  2u| 
directly  -  with  the  preferred  values  of  ka  and 


f  2 


sin  4u  +  16f  (7  -  9f)  sin  2u) 


9 


k. 

l 


076) 


as 


(177) 


4.9  Perturbation  in  U 

From  (20  and  (3)  we  obtain  at  once,  on  combining  (149)  with  (174)  and  (150) 
with  (175),  the  first-order  secular  result  for  Uj  ,  as  given  by  (111),  and  the 
periodic  term  given  by 


39 


Uj  -  -  A  [6(2  -  5f )  sin  2U  -  5  [6(26  -  33f)  sin  v 

+  3(4  -  Ilf)  sin  (u  +  5)  -  7(4  -  Ilf)  sin  (2u  +  v)]|  ,  (178) 

which  is  equivalent  to  (74) . 

The  expression  for  U2  is  given  by  combining  (154)  and  (176);  thus 

U2  =  *V  |(6  -  f  -  1 3f 2)  sin  4u  -  2  [f(40  -  33f)  -  12ka(2  -  5f) 

-  30kif(l  -  f)]  sin  2u}  .  (179) 

With  the  preferred  values  of  k  and  k.  , 

3  X 

U2  =  ^  |(6  -  f  -  1 3f 2)  sin  45+2  (24  -  106f  +  93f2)  sin  2u|  .  (180) 

Again,  the  general  expression  for  U2  is  given  by  combining  (153)  with  the 
second-order  component  of  An  ,  given  by  (174);  thus 

D2  =  [(432  -  1 052f  +  637f2)  -  24(3  -  4f )  (4kfl  +  3kR)  -  192k£f  ( 1  -  f)  +  48k2J 

.  (181) 

When  U.  is  zero,  through  k  having  its  preferred  value,  it  is  obviously 
1  •  n 

desirable  for  U2  to  be  zero  as  well,  and  (181)  indicates  that  this  requires 

k2n  =  "  &  {fO00  -  131  f )  -  96ka(3  -  4f)  -  192k.f(l  -  f)|  .  (182) 

With  the  preferred  values  of  ka  and  k^  ,  this  gives 

k,  -  A  (288  -  724f  +  515f2)  .  (183) 

2n 

The  associated  value  of  p2  is  given  by  (168).  Taking  k2n  as  given  by 
(182)  we  find 

P2  =  ^  | (324  -  8 1 4f  +  533f2)  -  120ka  (3  -  4f)  +  72k2  -  96kjf  (1  -  f)  -  36k2aj 

.  (184) 

but  with  the  preferred  values  of  ka,  k^  and  k2a  we  get  a  tremendous  simplifi¬ 
cation,  viz  to 

u2  =  A  f  (4  -  19f)  .  (185) 

It  may,  however,  be  considered  undesirable  to  tamper  with  the  Kepler  rela- 
tion,  in  which  case  y2  must  be  zero  and  we  cannot  force  U2  to  be  zero.  But, 
for  (168)  to  give  zero,  we  require 
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k0  =  rV  (5£2  +  16k2  -  24k  k  +  18k2  -  8k„  )  (186) 

2n  a  a  n  n  2a 

and  then  (181)  gives 

02  =  An  |(2 1 6  -  526f  +  34  If2)  -  12(3  -  4f )  (4kg  +  3kn)  -  96k. f()  -  f) 

+  72k2  -  108k  k  +  81k2  -  36k„  >  .  (187) 

a  a  u  n  2a  J 

Instead  of  substituting  our  preferred  value?  of  k  and  k  in  (’86),  it 

'  an 

is  instructive  to  set  k  =  k  =  jh  ,  with  k  given  by  (117),  ie  appropriate 

an  2a 

to  a  =  a'  ,  since  we  c<*u  then  interpret  il  a-  the  exact  constant  n'  .  li«. 
resulting  expression  for  k2n  is  given  ay 

k2n  =  "  ^  ^(4°  "  U8~l  “  !5^2)  “  72ki^  (1  “  b]  ;  0  88) 

with  the  preferred  k.  , 

l 

k„  =  -  ts  (40  -  1 20f  +  57f2)  .  (189) 

2n 

5  COMPARISONS  WITH  OTHER  AUTHORS'  RESULTS 

Of  the  various  papers  mentioned  in  section  ’ ,  comparisons  will  be  made 

only  for  Refs  22,  25  and  26.  The  two  papers  ftom  the  French  school,  one  by 

25  22 

Berger  and  Walch  and  the  other  by  Bretagnon  ,  are  considered  together  in 

section  5.2.  Since  they  both  give  results  for  e  and  w  ,  rather  than  £  and 

n  ,  suitable  formulae  for  comparison  are  first  derived  in  section  5.1.  The 

26 

remaining  paper,  by  Kinoshita  ,  is  considered  in  section  5.4  after  section  5.3 
has  given  results  needed  for  the  comparison,  namely,  for  the  quantities  fa.  and 

r~7  27 

•'av..  -  e  )  cos  l  . 

5 . 1  Second-order  analysis  for  e  and  m 

Formulae  for  second-order  perturbations  in  the  elements  e  and  u>  are 
necessarily  much  more  complicated  than  those  in  £  and  n  since  we  can  no 
longer  confine  ourselves  to  the  lowest  power  of  e  in  the  expansions.  Wa  use 
our  results  for  £  and  n  to  derive  the  formulae  we  require  for  e  and  u>  , 
employing  the  notation  ej  ,  e^  etc  as  in  (112)  and  (113). 

The  first-order  formulae  have  already  been  given,  viz  by  (77) ,  (78)  and 
(79)  for  ej,  oij  and  Uj  respectively.  It  is  convenient  to  split  ej  and  Wj 
according  to  the  formulae 


=  el0  +  e  en 


(190) 


4! 


and 

<*>,  =  e  1  (a>10  +  e  uijj)  ,  (191) 

where  e^,  e^,  and  to ^  j  are  all  independent  of  e  .  Since  we  shall 

require  two  components  for  both  e2  and  ,  it  is  convenient  to  split  these 
quantities  in  a  similar  way,  writing 

e2  =  5  1  (e20  +  *  e21)  0^2) 

and 

-2 

“2s5  ^“20  +  5  w21^  *  (193) 

It  will  be  observed  that  the  singularity  underlying  the  use  of  e  and  w  has 

made  e„  of  order  e  1  ,  something  not  intuitively  to  be  expected  perhaps,  while 

.  ^  _-2 
making  of  order  e 

We  shall  see,  in  (205)  and  (2 1 0)  following,  that  e^Q  and  Ujq  in  fact 
say  nothing  about  true  second-order  perturbations,  since  they  are  constrained  to 
algebraic  values  related  only  to  ej  and  Wj  .  All  the  information  is  in  e2j 
and  ,  which  is  why  two  sets  of  terms  are  required  for  each  of  e 2  and  u>2  • 

We  start  by  substituting  in  (1)  and  (2)  the  expansions 

e  =  e  +  K  ej  +  K2  e2  (194) 

and 

u>  =  S3  +  K  Wj  +  K2  e2  ,  (195) 

and  relating  the  resulting  expressions  to  the  corresponding  expansions  of  £ 
and  n  .  This  gives 


= 

e  1  cos 

3  - 

e  Wj  sin  £>  , 

096) 

nl  = 

e j  sin 

3  + 

e  u)j  cos  u>  , 

0  97) 

?2  = 

(e2  "  i 

I  0)^ 

cos  Si  -  (e  u2 

+  e!  V 

sin  Si 

(198) 

and 

n2  = 

(e2  -  * 

e  Wj 

')  sin  u  +  (e  u2 

+  el  “l> 

cos  S3 

(199) 

The  formulae 

for  ?2 

and  0^2 

are  of  course  known,  whence 

E2  and 

02  are  also 

known,  where 

<N 

w 

= 

£2  cos  u  +  n2  sin  u 

(200) 

and 

°2 

a 

-  ?2  sin  “  +  ^ 

cos  u 

(201) 
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The  definitions  involved  in  (200)  and  (201)  will  be  extremely  useful. 
The  basis  for  the  notation  is  that  the  right-hand  sides,  when  the  2-suffices 
are  omitted,  give  e  and  zero  respectively.  It  is  now  immediate  that 


and 


E2  =  e2  -  Je  <o2  (202) 

02  =  i  u>2  +  e}  Uj  .  (203) 


We  can  now  substitute  (192)  and  (191)  in  (202)  to  obtain 

e2  =  e  1  (e2Q  +  e  e2J)  -  Je  1  (oo j Q  +  e  u)n)2  .  (204) 

Since  there  are  no  e  1  terms  in  E2  -  for  economy  of  space  its  expression  is 
not  given  -  it  follows  from  (204)  that 

e20  "  <205> 

and  this  is  the  first  'constraint*  referred  to  earlier.  An  explicit  formula 
follows  at  once  from  (79),  from  which  the  and  ^  components  of  are 

defined  by  (191).  We  get 


e 


20 


6(24  -  72f  +  6 1 f 2 )  cos  2v  -  240fh  cos  2u  +  168fh  cos  2(u  + 

-2  -2  -2 
+  9f  cos  2(u  +  u>)  -  42f  cos  4u  +  49f  cos  2(2u  +  v) 

+  72fh  cos  2«  -  2(72  -  216f  +  191f2) j  . 


v) 


(206) 


From  identification  of  the  e°  terms  of  (204)  we  'jee  that  the  other 
component  of  e9  is  given  by 


'21 


=  E2  +  “10  “ll 


(207) 


On  evaluation  of  the  right-hand  side,  we  get 

e21  =  ^  {9  C144  "  336^  +  2,,f2  ”  128kah  “  96kif(l  -  f)]  cos  v 

-  9(16  -  48f  +  39f2)  cos  3v  +  6f  [l34  -  189f  -  48k&  +  24ki ( 1  -  f)]  x 
x  cos  (u  +  55)  +  8f  [28  -  65f  -  84kfl  +  42^(1  -  f)]  cos  (2u  +  v) 

-  192fh  cos  (2u  +  3v)  -  36fh  cos  (v  -  2ui)  +  36f2  cos  (3u  +  u) 

-  81 f 2  cos  (4u  +  v)  —  63f 2  cos  (4u  +  3v)  +  6k^  [l2h  sin  v 

-  3f  sin  (u  +  u)  +  7f  sin  (2u  +  v)T  1  .  (208) 
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But  the  preferred  k^  -  see  (93)  -  is  a  multiple  of  sin  22  ,  causing 
terms  to  combine  with  earlier  terms  of  (208) . 

In  the  same  way  we  can  substitute  (190),  (191),  and  (193)  in  (203)  to 

obtain 


°2  "  5’'  (u)20  +  5  “215  +  S~'  (e10  +  5  ell)(<°10  +  5  “ll}  * 

_-l 


(209) 


from  which,  since  02  contains  no  e  terms. 


di 


20 


e10  “l0 


(210) 


This  is  the  other  'constraint'  and  it  leads  to  the  explicit  formula 


20 


=  -  f88  |b(24  -  72f  +  6 1 f ^ )  sin  2v  +  168fh  sin  2(0  +  v) 


-  9f2  sin  2(0  +  2)  +  49f2  sin  2(20  +  v)  -  72fh  sin  22^ 


(211) 


Also,  from  identification  of  the  5°  terms  of  (209)  we  have 


21 


“  °2  ”  e10  “ll  ell  “lO  » 


(212) 


evaluation  of  which  gives 

u>2,  -  ^6  [96  -  216f  +  101  f 2  -  96kah  -  72k.f(l  -  f)]  sin  v 

-  9(16  -  48f  +  39f2)  sin  3v  -  18?  [lO  -  13f  -  8ka  +  4^(1  -  f)]  x 
x  sin  (u  +  2)  -  2f  [70  -  59f  +  168ka  -  84ki(l  -  f)]  sin  (20  +  v) 

-  192fh  sin  (2u  +  3V)  -  36fh  sin  (?  -  22)  +  27f2  sin  (30  +  GJ) 

-  I59f2  sin  (40  +  7)  -  63f2  sin  (40  +  3v)  -  3kg  [l2h  sin  v 

-  3f  sin  (u  +  2)  +  7f  sin  (2u  +  v)T  -  3k  [l2h  cos  v 

+  3f  cos  (u  +  2)  +  7f  cos  (20  +  v) 


]}. 

(213) 


Again  the  ke  and  k^  terms  combine  with  earlier  terms,  if  preferred  values  are 
introduced . 

5.2  Comparison  with  Berger  and  Walch,  and  with  Bretagnon 

Refs  25  and  22  both  tabulate  large  numbers  of  terms  in  a  compact  manner,  of 
which  only  a  small  number  are  involved  in  the  comparisons.  A  convenient  reference 
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system  is  desirable,  and  the  following  is  adopted:  BW  P.l  refers  to  line  1  of 
page  P  of  Ref  25,  while  B  P.pqn  refers  to  the  particular  line  of  the  appro¬ 
priate  table  on  page  P  of  Ref  22  that  is  designated  by  the  given  values  of 
p,  q  and  n  .  Both  papers  in  fact  tabulate  numerical  coefficients  against  the 
trigonometric  arguments  (poo  +  qM)  and  cover  (in  a  single  line  of  tabulation) 

more  than  just  the  lowest  power  of  e  that  is  relevant  here  (the  powers  of  e 

2. 

increase  in  steps  of  two) .  The  index  n  relates  to  powers  of  f  (=  sin  i)  that 
are  allocated  different  lines  of  Ref  22  though  appearing  on  the  same  line  of 
Ref  25. 

To  get  agreement,  it  will  be  found  that  we  require  k  =  k.  =  0  at  all 

3  X 

times  in  Ref  25,  but  various  values,  that  will  be  stated,  for  Ref  22.  The 
comparison  proceeds  through  the  elements  in  the  usual  order,  except  that  e  and 
a)  are  held  over  till  last. 

Starting  with  a  ,  then,  the  second-order  formula  given  here  as  (121) 

(with  =  0  )  is  in  agreement  with  BW  122.18  (for  cos  4u)  and  (with 

k  =  k.  =  0  )  BW  122.15  (for  cos  2u) .  It  also  agrees  with  B  16.444  , 

B  16.222  and  B  16.224  ,  so  long  as  k  is  set  to  -2h  and  k.  to  1 . 

3  1 

For  i  ,  the  second-order  formula  given  here  by  (127)  (with  k,^  -  0  )  is 

in  agreement  with  BW  122.47  and  BW  122.45  .  It  also  agrees  with  B  22.440  , 

B  22.442  ,  B  22.220  and  B  22.222  ,  so  long  as  k  is  set  to  -2h  and  k. 
to  -2  . 

For  £2  ,  we  have  a  secular  effect  as  well  as  short-periodic  terms.  The 

contribution  to  0  ,  represented  here  by  (133),  is  the  same  as  is  given  by 

BW  110.20  ,  so  long  as  we  set  k  (as  well  as  k  and  k .  )  to  zero;  the  same 

n  3  i 

value  is  indicated  by  Ref  22  through  reference  to  Brouwer's  paper  (see  the  second 
equation  on  page  394  of  Ref  11).  For  the  periodic  effect,  the  formula  given 
here  as  (134)  is  in  agreement  with  BW  122.53  and  BW  122.51  ;  for  agreement 
with  B  22.440  ,  B  22.442  ,  B  22.220  and  B  22.222  we  must  set  ka  to  -ljh 
and  k^  to  -1  . 

For  U  ,  we  have  a  secular  effect  that  (for  an  unmodified  Kepler  relation) 
is  given  here  by  (187),  and  this  agrees  with  the  combination  of  BW  110.25  and 
BW  110.30  if  k^  and  k^a  are  set  to  zero.  Bretagnon  (page  II  of  Ref  22) 
gives  the  BW  110.30  formula,  ie  for  >  explicitly,  on  the  grounds  that 
there  was  an  error  in  Brouwer's  paper  (last  equation  on  page  393  of  Ref  P);  as 
we  shall  see  in  sections  5.4  and  9.7,  the  discrepancy  amounts  to  an  implicit  use 
of  a  non-zero  k2a  by  Brouwer  -  in  giving  the  explicit  formula,  Bretagnon 
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unfortunately  introduces  two  real  errors,  since  he  writes  J2  instead  of 
and,  less  obviously,  halves  the  main  term,  which  should  contain  27/4  and  not  27/8 
For  the  periodic  effect,  the  formula  given  here  as  (179)  is  in  agreement  with 
the  combination  of  BW  123.14  with  BW  123.39  and  of  BW  123.2  with  BW  123.27 
again,  the  three  combinations  given  by  B  28.440,  B  28.442  and  B  28.444  are 
also  easily  verified,  and  so  are  the  three  given  by  B  28.220,  B  28.222  and 
B  28.224  ,  so  long  as  we  set  kfl  to  -l^h  and  k^  to  . 

We  complete  the  comparisons  by  looking  at  e  and  u>  together,  calling  on 

results  from  section  5.1.  We  have  results  of  two  orders  of  e  to  check,  as  has 
been  explained.  The  lower-order  formulae  are  given  here  by  (206)  and  (211),  and 
they  are  in  immediate  agreement  with  terms  given  by  BW  122.22,  BW  122.27, 

BW  122.29,  BW  122.35,  BW  122.37  and  BW  122.39  ,  for  e2Q  ,  and  BW  122.56, 

BW  123.4,  BW  123.12  and  BW  123.16  ,  for  >  constant  terms  being  ignored; 

the  same  results  are  given  by  B  16.  20,B16.  22,  B  16.  24,  B  16.222,  B  16.224, 

B  16.242,  B  16.244,  B  16.424,  B  16.444  and  B  16.464  ,  for  e2Q  ,  and  B  28.  20, 

B  28.  22,  B  28.  24,  B  28.242,  B  28.244,  B  28.424  and  B  28.464  ,  for  w2Q  . 

The  higher-order  formulae  are  given  here  by  (208)  and  (213),  but  two 

corrections  are  needed  before  they  give  agreement  with  BW  122.21,  BW  122.23, 

BW  122.26,  BW  122.28,  BW  122.30,  BW  122.33,  BW  122.36,  BW  123.38  and  BW  123.40 

for  e2J  ,  and  BW  122.55,  BW  122.57,  BW  123.1,  BW  123.3,  BW  123.5,  BW  123.9, 

BW  123.13,  BW  123.15  and  BW  123.17  ,  for  w2]  .  First,  allowance  has  to  be 

made  for  carry-over  terms  from  the  lower-order  formulae;  these  arise  because  the 

argument  is  u  here  but  U  in  Ref  25.  Secondly,  non-zero  values  have  to  be 

set  for  kg  and  k^  ;  we  get  the  desired  agreement  if  we  take*  2f  cos  2w  +  8fi 

for  k  and  -2f  sin  2 u  for  k  .  Ref  22  gives  the  same  results  (assuming  the 

e  u) 

same  kg  and  k^  ),  the  e2]  terms  being  from  B  17.  10,  B  17.  12,  B  17.  14, 

B  17.  30,  B  17.  32,  B  17.  34,  B  17.212,  B  17.214,  B  17.2-12,  B  17.2-14,  B  17.232, 

B  17.234,  B  17.252,  B  17.254,  B  17.434,  B  17.454  and  B  17.474  ,  but  we  have  to 

take  k  =  -Jh  and  k.  =  0  when  pq  =  1  ,  k  =  - 1  Ih  and  k.  =  2  when  pq  =  21 
and  k^  =  -£h  and  k^  =  $  when  pq  =  23  .  For  the  u)^  terms  agreement  is 
obtained,  similarly,  from  B  29.  10,  B  29.  12,  B  29.  14,  B  29.  30,  B  29.  32, 

B  29.  34,  B  29.212,  B  29.214,  B  29.2-12,  B  29.2-14,  B  29.232,  B  29.234,  B  29.252, 

B  29.254,  B  2°. 434,  B  29.454  and  B  29.474  ,  subject  to  the  same  assumptions 

concerning  k^  and  k^  . 

*  These  va’ues  of  ke  and  k^  are  actually  0(e)  truncations  of  the  ke  and 
kw  that  are  implicit  in  Refs  22,  25  and  26;  they  are  defined  such  that  the 
time-averages  of  6e  and  6u  are  zero,  and  the  full  expressions  can  be 
obtained  from  section  9. 
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This  paper  has  so  far  not  presented  a  formula  for  the  second-order  secular 
perturbation  in  10  ,  since  this  perturbation  is  O(Ke)  in  its  effect  on  satellite, 
position  (apart  from  the  effect  in  combination  with  M  ,  which  is  allowed  for  by 
0  ),  as  explained  in  section  3;  5L  is  0(K)  in  fact,  as  we  shall  see,  whereas 

-2  1  -j 

W2Q  contributes  an  0(K  e  )  term  to  6io  and  contributes  an  0(K  e  ) 

tern.  However,  it  is  perhaps  worth  deriving,  and  checking,  a  formula  for  ai2  , 
so  that  the  secular  perturbations  are  all  known  to  0(Ke)  -  see  also  section  9.6. 

We  can  start  with  Merson's  formula  for  “sec  >  given  by  equation  (A-27)  of 
Ref  7.  This  yields 

S2  =  H  nf  (4  -  5f )  ;  (214) 

the  underlying  values  of  the  k-constants  here  are  the  Kozai  values,  viz 

k  =  h,  k.  =  0  and  k  ■  lh  .  From  consideration  of  the  (first-order) 

3  x  n 

expression  for  to  ,  given  by  (78) ,  we  can  obtain  the  general  contribution  of 
the  k's,  from  which  the  general  formula  must  be 

3.  =  ^n  ((288  -  676f  +  395f2)  -  1 2 (4k  +  3k  )  (4  -  5f)  -  120k. f(l  -  f)l  . 

/I  an  l  J 

.  (215) 

We  now  put  k&  »  k^  =  kR  =  0  ,  to  obtain  a  result  which  agrees  with  BW  110.25. 

In  terms  of  the  preferred  constants  of  the  present  paper  we  derive 


m2  =  -  **gn  (192  -  476f  +  325f2)  .  (216) 


If  'fz  is  denoted  by  a  ,  the  (second-order)  formula  for  c*2  is  easily 
obtained. 

The  basic  identity  is 

(5  +  ic  0]  +  K2  a2)2  =  5  +  ic  at  +  ic2  a2  , 

which  yields  the  first-order  relation 
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Taking  and  a 2  from  (65)  and  (121),  we  get 

a  =  **ga  <5f^  cos  4u  +  4f  ^2(5  -  9f)  -  9k  +  6k  (1  -  f)]  cos  2u 

-  «f2  ♦  ^  -  8k2a>}  ' 

.  (217) 

2  2. 

We  also  need  to  know  how  the  quantity  a(l  -  e  )  cos  i  varies,  or  rather 

2 

its  square  root  which  we  denote  by  6  .  Now  a(l  -  e  )  is  simply  p  ,  the 
semi-latus  rectum  of  the  osculating  elliptic  orbit,  and  it  is  an  elementary 
result  that  pp  =  h  ,  where  h  is  the  angular  momentum  of  the  satellite's 
motion  relative  to  the  centre  of  the  earth.  Clearly  h  acts  normally  to  the 
orbital  plane,  so  h  cos  i  is  the  angular  momentum  about  the  earth's  axis.  But 
this  is  constant  for  the  force  associated  with  any  axi-sytmetrio  potential,  ie 
for  any  zonal  harmonic.  Also  h  cos  i  =  3  ,  bo  our  result  is  that  3  ,  like  a' 
and  n'  ,  is  an  absolute  constant  of  the  motion  in  the  J2  field. 

5.4  Comparison  with  Kinoshita 

Appendix  B.l  of  Ref  26  lists  the  second-order  short-period  perturbations 
2 

due  to  •  There  are  many  terms  in  each  of  the  five  formulae  given,  and  a 
convenient  reference  system  is  to  let  KB-P.t  denote  term  t  of  Kinoshita 's 
page  B-P  . 

The  five  formulae  relate,  in  order,  to  elements  that  identify  with  £,  n, 

U,  fl  and  i/pa  in  the  notation  of  the  present  paper.  A  sixth  formula  is 
unnecessary  since  Kinoshita's  sixth  element  is  y|pa(l  -  e^)^  cos  i  »  an<*  we 
have  seen  that  this  is  unperturbed.  To  get  agreement  we  again  require 
kfl  =  k^  =  0  ,  with  also  zero  for  secular  quantities.  (Like  Berger  and 

Walch,  and  Bretagnon,  Kinoshita  expresses  results  in  terms  of  M  ,  rather  than 
v  ,  so  that  his  mean  elements  are  true  time  averages.) 

We  start  with  S2  »  then,  finding  that  the  second-order  formula  given 

here  as  (141)  is  in  agreement  with  KB-2.8,  KB-3.4  and  KB-4.1  .  For  n2  * 

formula  (146)  here  agrees  with  KB-5.6,  KB-6.3  and  KB-6.10  .  For  U2  ,  (179) 

here  agrees  with  KB-8,7  and  KB-9.1  ;  also  U2  given  by  (187)  agrees  with 

the  combination  of  the  fourth  and  fifth  formulae  of  KB-53  ,  so  long  as  we  set 

k_  to  (12  -  30f  +  25f^)/I2  ,  a  value  that  can  be  obtained  empirically  but 

2a 

which  is  confirmed  by  the  formula  for  v'ya  .  (The  fifth  formula  of  KB-53  ,  on 
its  own,  agrees  with  S2  »  38  given  by  (215)  with  every  k  set  to  zero.)  For 
,  (134)  here  agrees  with  KB-I0.9  and  KB-10.14  ;  also  given  by  (133) 
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agrees  with  the  sixth  formula  of  KB-53  .  For  v'ya  ,  (217)  here  agrees  with 
KB-12.8,  KB-12.13  and  KB-13.2  ,  if  is  assigned  precisely  the  'empirical' 

value  just  quoted.  (Translated  into  k2n  ,  by  multiplication  by  -1$  in  accord¬ 
ance  with  ( <  86) ,  this  k2a  fully  resolves  the  Brouwer-Bretagnon  discord 
remarked  in  section  5.2.) 

This  completes  the  tie-up  between  the  present  results  and  those  of  preced¬ 
ing  papers,  apart  from  some  results  for  (untruncated)  secular  and  long-periodic 
perturbations  which  are  derived  and  verified  in  section  9. 

6  SECOND-ORDER  PERTURBATIONS  IN  POSITION 

To  obtain  perturbations  in  the  cylindrical  polar  coordinates  described  in 
section  2.5,  it  is  essentially  a  matter  of  combining  existing  results,  no 
further  integration  being  required.  As  a  simplification,  we  throughout  set 
k  ,  k  ,  k  and  k.  to  their  preferred  values,  as  given  by  (89),  (90),  (91) 

3  ^  I]  1 

and  (99),  since  this  preference  was  based  on  the  resulting  simplicity  in  the 
first-order  formulae  for  perturbations  in  coordinates. 

It  is  convenient  to  give,  in  advance,  some  of  the  formulae  that  will  be 
required.  Thus,  neglecting  0(,5)  terms  in  each  case, 


*1 

cos 

u 

+  nl 

sin 

5  * 

h  +  &f  cos 

ho 

Cl 

(218) 

sin 

u 

-  "l 

cos 

u  = 

-  if  sin  25 

9 

(219) 

e2 

cos 

u 

+  n2 

sin 

u  = 

-  &f  i3(l 

-  f)  cos  40  -  2(10  - 

Ilf)  cos  25  -  3(1  - 

7»} 

\ 

(220) 

sin 

u 

"  n2 

COS 

u  = 

rh  |3f(io 

-  1 3f  ) 

sin 

45  +  2(36 

-  40f  -  f2)  sin  2s| 

(221) 

and 

0? 

-  r 

2^ 

V 

)  sin 

25  - 

■  25jfi 

1  cos  25  = 

-  iVf 

(5f 

sin  4u  + 

12h  sin  25) 

(222) 

These  formulae  follow  at  once  from  (70),  (72),  (142)  and  (147). 

6 . 1  Perturbation  in  r 

We  require  formulae  for  rj  and  r 2  >  such  that 

r  =  ?  +  K  rj  +  K2  r2  ,  (223) 

where  r  is  derived  from  mean  elements  by  application  of  the  standard  algorithm 
of  section  2.5,  whilst  rj  and  r2  are  purely  periodic. 


100 
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From  (84),  it  follows  that 

r  -  | a  +  Kaj+  K2a2|jl  -  (l  +  K£,  +  K2£2)(cos  U  -  KUj  sin  U) 

-  (n  +  Krij  +  K2n2)(sin  U  +  KUj  cos  U)  +  [(£  +  KC , )  sin  U 

"  (n  +  Krij)  cos  u]2}  ,  (224) 

where  only  such  terms  are  retained  as  will  actually  be  needed. 

On  expanding  (224)  and  comparing  with  (223),  we  obtain  the  first-order 
formula 

Tj  -  (r/3)  a!  -  a  [(£]  cos  U  +  sin  U)  -  e  sin  M  (Uj  +  2£j  sin  U  -  2n ^  cos  U)]  , 

.  (225) 

which  reduces  to  (85),  and  hence  (88),  in  view  of  (53).  The  formula  for  r2  , 
similarly,  is 

r2  =  (r/a)  a2  -  a]  (£j  cos  0  +  n j  sin  U)  -  a  [(£,>  cos  U  +  n2  sin  U) 

-  (£j  sin  U  -  n,  cos  U)(£j  sin  U  -  Pj  cos  U  +  Uj)]  .  (226) 

We  can  replace  U  by  u  of  course,  since  we  are  ignoring  0(Ke)  pertur¬ 
bations.  Then  r2  is  a  combination  of  terms  given  by  (122),  (65),  (218),  (220), 
(219)  and  (74),  the  result  on  reduction  being 

r2  =  -  7^a  |f2  cos  40  +  2f  (26  -  31  f)  cos  20  +  72h  -  f2  -  72k2a}  . 

.  (227) 

Choice  of  k2a  ,  as  foreshadowed  by  (123),  now  gives  our  final  result,  viz 

r2  ~  ~  ^Vaf  | f  cos  4u  +  2  (26  -  31f)  cos  2u  |  .  (228) 

6.2  Perturbation  in  u 

We  require  formulae  for  Uj  and  u2  ,  such  that 

u  =  u+Kuj+K2u2  ,  (229) 

where  Uj  and  u2  are  purely  periodic. 
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From  (94)  it  follows  that 

u  =  (U  +  KU,  +  K2U2)  +2(1  +  KC,  +  K2C2)(sin  U  +  KU,  cos  U) 

-  2(n  +  Kn ,  +  K2n2)(cos  0  -  KU I  sin  U)  +  1 1  {[(I  +  K£,)2 

-  (n  +  Kn,)2]  Sin  2U  -  2(1  +  KC,)(n  +  Kn,)  cos  2u}  ,  (230) 

with  only  the  needed  terms  retained. 

On  expanding  (230)  and  comparing  with  (229),  we  obtain  the  first-order 
formula 

Uj  *  U,  +  2  (5,  sin  U  -  n,  cos  U)  +  $e  [4Uj  cos  M  +  5£,  sin  (Q  +  v) 

-  5n,  cos  (u  +  v)]  ,  (231) 

which  reduces  to  (95)  in  view  of  (53).  The  formula  for  u2  ,  similarly,  is 

u2  *»  U2  +  2  (C2  sin  U  -  n2  cos  U)  +  2U j  (C,  cos  U  +  n,  sin  U) 

+  li  [(cj  -  n2)  Sin  2U  -  2€,n,  cos  2uJ  .  (232) 

We  replace  U  by  u  and  obtain  u2  as  a  combination  of  terms  given  by 

(180),  (221),  (74),  (218)  and  (222),  the  result  being 


u2  »  rh  {(,8  "  3f  -  17f2)  sin  4G  +  2  (72  -  170f  +  97f2)  sin  2g|  .  (233) 

6.3  Perturbations  in  r',  u*  and  c 

The  first-order  results,  given  in  section  3,  were  based  on  formulae  (43)  and 
(44),  which  are  not  valid  when  we  proceed  to  second  order.  Thus  we  require  a 
fresh  start,  based  on  the  formula  (45)  that  effectively  defines  the  cylindrical 
polar  coordinates  r',  u'  and  c  . 

T 

Now  (x  y  z)  is  given  by  both  (41)  and  (45),  so  an  exact  equation 


for  r',  u* 

and 

c  is 

(r'  cos  u' 

r' 

sin  u* 

■  R, (i)  R3(0)  R3(-0)  Rj (— i)  (r  cos  u  r  sin  u  o)T  .  (234) 

We  want  to  use  (234)  as  source  for  perturbation  formulae  that  are  correct  to 
second  order.  Now  R3(S)  R3(-ft)  ■  Rj(Q  -  Q)  and  of  course  ft  -  5  *  5ft  by 


l 


i 

j 
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the  basic  definition  (33).  We  write  5  in  place  of  ,  for  convenience,  and 
also  write 


1  =  i  -  i  (  =  6i)  , 

with  sine  and  cosine  of  i  and  i  denoted  by  a,  a,  s  and  a  .  Then  the 
second-order  approximation  of  (234)  is 


r  cos  u 


1  -  $6 


1  -  Ji2  -  \c2 &2 
2 

t  +  $sc?6 


-t  +  $806 

2  2  2 

1  -  $i  -  Js  SI 


r  cos  u 


r  sin  u 


The  individual  equations  may  be  written  as 


r*  cos  u'  =  r  (cos  u  -  06  sin  u  -  $6  cos  u)  , 


r  sin  u 


r  ^sin  u  +  56  cos  u  -  $  2  +  c2 &2~\  sin  u) 


(235) 


(236) 


c  =  r(i  sin  u  -  s6  cos  u  +  Jsc6  sin  u) 


(237) 


To  first  order  we  again  see  that  r'  =  r  and  recover  the  results  given  by 
(102)  and  (100)  of  section  3. 

To  second  order,  we  still  have  identity  between  r'  and  r  ,  consequent 

on  our  choice  of  k.  ,  as  remarked  at  the  end  of  section  2.5.  Thus  we  need  not 
i 

consider  r'  as  a  separate  quantity,  but  note  that  we  would  have  to  do  so  if 
.  2  2 

we  were  concerned  with  0(Ke  )  or  0(K  e)  perturbations. 

We  can  now  obtain  the  second-order  perturbation  in  u'  ,  using  either 
(235)  or  (236)  -  as  a  check  we  use  both .  From  (235)  we  have 


cos  u'  =  cos  u  -  (3  -  3t)  6  sin  u  -  Ji  cos  u  , 


(238) 


where  6  =  KT2  j  +  K  and  i  =  Kij  ,  while  from  (236)  we  have 

2  2  2  . 

sin  u’  =  sin  u  +  56  cos  u-$(t  +  a  6  )  sin  u 


(239) 


But  if  u  denotes  u'  -  u  ,  we  have  the  Taylor  expansions 


cos  u'  =  cos  u  -  u  sin  u  -  $u  cos  u 


.  ,  .  ,  i. 

sin  u  =  sin  u  +  u  cos  u  -  ju  sin  u 
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We  identify  the  two  expressions  for  cos  u'  ,  and  likewise  for  sin  u'  , 
expanding  u  according  to  the  formula 

u  =  K  u  +  K2  u2  ,  (240) 

where  Uj  and  u2  are  required.  It  is  immediate  that 

u,  -5a,  (241) 

but  we  have  two  formulae  for  u2  ,  viz 

2  2 

(u2  -  efi2)  sin  u  =  -  si , ,  sin  u  +  £(ft]  ~  Uj)  cos  u  (242) 

from  the  cos  u'  identity,  and 

-  2  2  2  2 

(u2  -  ofi2)  cos  u  =  J  (u,  -  i,  -  c  flj)  sin  u  (243) 

from  the  sin  u'  identity.  On  substitution  for  ,  i^  and  ,  we  get  the 
same  formula  in  both  cases,  viz 

U2  "  3fi2  =  ~  ^  ( 1  ~  f)(sin  45+2  sin  25)  .  (244) 

From  (136),  the  formula  for  fi2  ,  we  now  get 

u2  =  -  ^  (1  -  f) |(6  +  5f )  sin  45  +  2(24  -  39f)  sin  25  J-  .  (245) 

For  u2  ,  it  remains  to  combine  (233)  and  (245),  the  result  being 

u2  =  -  T1!?  jf  sin  45  -  (19  -  20f)  sin  2u|  .  (246) 

We  finally  obtain  Cj  and  c2  ,  where 

c  (=  6c)  =  K  c,  +  K2  c2  ,  (247) 

from  (237).  We  consider  this  equation  in  the  form 

-  -2-2 
c/r  -  K  ( i ,  sin  u  -  e  fij  cos  u)  +  K  (i2  sin  u  -  s  &2  cos  u  +  ,  sin  u) 

.  (248) 

Now  from  (66)  and  (68) ,  identifying  u  and  u  only  in  the  e  term, 

i,  sin  u  -  ft,  sin  i  cos  u  -  j  sin  i  cos  i  |[sin  u  -  sin  (2u  -  u)] 

+  l^e  [2  sin  (0  +  v)  -  3  sin  u)]  j 

.  (249) 


1 
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Thus  we  have  recovered  (100),  except  that  there  is  a  ’second-order  carry-over' 
given  by  the  first  term  of  (249),  which  can  be  written  as  J (u  -  u)  sin  2i  cos  u  . 
The  second-order  terms  now  combine  to  give 

_  2  — 

c^/r  =  i2  sin  u  -  sin  i  cos  u  +  Jftj  sin  2i  sin  u  +  Juj  sin  2i  cos  u  , 

.  (250) 

which  from  (128),  (135),  (68)  and  (97),  replacing  u  by  u  ,  reduces  to 

c2  =  -  fa  r  sin  2i  ^8f  sin  3u  -  (23  -  49f  -  k2.)  sin  u|  .  (251) 

The  reason  for  choosing  k^  according  to  (129)  should  now  be  apparent  and  the 
final  result  is  simply 

c2  =  -  T^irf  sin  2i  sin  3u  .  (252) 

It  is  worth  remarking  that,  as  with  c^  given  by  (100),  we  can  replace  the  factor 

r  by  either  r  or  a  if  desired.  It  no  longer  makes  sense  to  derive  an 

from  u2  ,  however,  by  multiplying  by  any  of  these  factors. 

6.4  Perturbations  in  instantaneously  fixed  RLC  directions 

The  basis  for  the  immediately  preceding  remark  (concerning  and  u2) 

is  that  a  natural  second-order  quantity  JZ-2  cannot  be  defined,  since  the  second 
coordinate  of  a  cylindrical  triple  is  unalterably  angular.  (This  point  has 
already  been  made  in  sections  1  and  2.5.)  Thus  our  fundamental  results  are  for 
r2>  u2  and  c2  ,  given  by  (228),  (246)  and  (252)  respectively. 

For  completeness,  however,  we  can  obtain  formulae  for  perturbations  in 
three  orthogonal  directions,  viz  the  RLC  directions  as  defined  in  section  2.5  by 
irstantaneous  values  of  the  mean  elements.  Denoting  coordinates  in  these  (fixed) 
directions  as  (x' ,  y’,  z')  ,  we  have 

(x' ,  y',  z')  =  £r'  cos  (u'  -  u) ,  r'  sin  (u'  -  u) ,  c]  ,  (253) 

where  (x',  y',  z')  =  (?,  0,  0) 

Introducing  the  formal  expansions  of  r' (  =  r  to  the  order  required)  and 
u' ,  with  corresponding  expansions  for  x',  y'  and  z',  we  find  at  once  that 

x]  =  r]  ,  yj  =  ruj  (of  (47))  and  zj  =  cf  , 

whilst  the  second-order  formulae  are 
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and 


r2  -  J  f  uj2  , 

(254) 

r  U2  +  r]  U1 

(255) 

=  c 


2  • 


(The  resemblance  of  (254)  and  (255)  to  familiar  formulae  of  dynamics,  immediately 
clear  on  replacing  r ^  by  Jr  eta ,  is  of  course  no  coincidence.) 

Evaluation,  based  on  (88),  (102),  (228)  and  (246),  gives 


X2  =  ~  s"taraf  [7f  cos  4u  +  16  (26  -  31  f)  cos  2u  +  fj 

and 

yl  =  -  risaf  [f  sin  4u  -  2  (19  -  20f)  sin  2u]  . 

7  NUMERICAL  CHECKS 

In  addition  to  the  formulae  checks  against  Refs  22,  25  and  26,  as  described 
in  section  5,  it  was  decided  to  make  two  numerical  checks  by  comparing  a  hand- 
computation  of  the  formulae  of  this  paper  with  the  results  of  a  numerical  inte¬ 
gration.  Both  these  numerical  checks  were  based  on  an  orbit  of  zero  eccentricity, 
63°  inclination  and  12h  period,  akin  to  that  of  NTS  2  (Ref  27),  with  a  greatly 
magnified  value  of  J2  and  all  other  sources  of  perturbation  suppressed. 

For  the  first  check,  J2  was  set  to  0.05  and  epoch  conditions  were  chosen 
to  make  mean  eccentricity  about  10  effectively  zero  for  most  purposes,  the 

.  .  ...  —  O 

corresponding  osaulattng  eccentricity  being  about  10  .  With  y  =  398602.0,  an 

ephemeras  was  generated  by  numerical  integration,  and  the  values  of  (osculating) 

a,  i,  U,  ?,  n>  U  noted  at  epoch  and  after  7  hours,  ie  210°  around  the  orbit. 

The  changes  in  the  elements,  and  also  in  r  and  u  ,  were  compared  with  values 

computed  from  the  formulae  of  sections  4  and  6.  The  result  of  this  check  was 

entirely  satisfactory,  the  discrepancy  for  each  comparison  being  of  order 
-3-9 

K  (  *=  81  x  10  for  the  magnified  )  as  expected.  In  particular,  with  our 

usual  definition  of  n  and  a  ;  such  that  h  is  given  by  U  (and  hence 

i  _2_3 

includes  u  )  and  k  is  given  by  (123),  we  have  n  a  =  398980.92  ,  in 

jLSL 

conformity  with  the  expression 

n2a3  =  yjl  +  AZ  [l 2 (6  -  7f)  +  Kf(4  -  19f)]j  (256) 

that  results  from  substitutions  of  (162)  and  (185)  in  (158). 
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The  second  check  was  more  elaborate  and  involved  use  of  the  NED  computer 

27  .  .  .  . 

program  in  a  specially  modified  form.  A  12-hour  ephemens,  of  interval  15  min, 

was  generated  by  numerical  integration,  on  this  occasion  with  J ^  set  to 
0.01082628  (ten  times  the  value  normally  used,  but  corresponding  to  a  value  of 
K  that  would  be  correct  for  a  close  earth  satellite),  p  =  398601 .3  and  osculating 
eccentricity  zero  at  epoch  (mid-point  of  the  12-hour  period),  the  corresponding 
mean  eccentricity  being  0.000254.  The  following  22  parameters  of  the  modified 
NED  were  fitted  to  the  29  points  (ie  87  coordinate  quantities)  of  the  ephemeris: 
a,  e,  i,  fl,  u,  0,  n,  H  and  14  coefficients  in  expressions  representing  perturba¬ 
tions  in  r,  u'  and  c  .  For  r  ,  the  perturbation  terms  involved  cos  2u  and 
cos  4u  ,  as  demanded  by  (88)  and  (228),  together  with  nt  cos  u  and  nt  sin  u  , 
which  are  required  because  NED  only  allows  for  w  as  a  component  of  n  (see  the 
remarks  on  suppression  of  ai  in  section  3).  For  u'  ,  the  perturbation  terms 
involved  sin  2u,  sin  4u,  sin  u  and  cos  u  ,  as  demanded  by  (102)  and  (246), 
together  with  nt  cos  u  and  fit  sin  u  .  For  c  ,  the  perturbation  terms 
involved  a  constant,  sin  2u,  cos  2u  and  sin  30  ,  as  demanded  by  (100)  and  (252). 

The  fitted  values  for  the  22  parameters  all  agreed  with  computed  values 
from  the  formulae  of  this  paper,  to  within  order  k\  **  0.82  x  10  ^  for  the 
magnified  J 2  ,  equivalent  to  22  mm  at  the  geocentric  distance  of  NTS  2).  The 
standard  deviation  for  the  87  coordinate  residuals  was  a  mere  0.7  mm. 

8  COMPLETION  OF  FIRST-ORDER  ANALYSIS  FOR  POSITION 

To  complete  the  analysis  of  section  3,  it  is  worth  giving  the  full 
(untruncated)  expressions  for  perturbations  in  the  cylindrical  polar  coordinates 
used  in  this  paper.  The  expressions  are  quite  compact  and  suggest  that  the 
correspondingly  complete  second-order  expressions  might  not  be  unduly  complicated. 

The  (accurate)  first-order  differential  expressions  for  the  cylindrical 

27 

coordinates  in  terms  of  elements  are 

6r  =  (r/a)  6a  -  a  6e  cos  v  +  (a/q)  e  6M  sin  v  ,  (257) 

-1  2.-1 

6u'  =  6a)  +  6fi  cos  i  +  q  5M  +  (2/q  )  (6e  sin  v  +  9  e  5M  cos  v) 

+  (e/2q  )  ^Se  sin  2v  +  q  e  6M  (3  +  cos  2v)J  (258) 

and 

6c  =  r  (6i  sin  u  -  612  sin  i  cos  u)  .  (259) 


In  these  expressions  we  substitute  for  6a  etc  from  equations  (55)  etc ,  except 
that  we  must  introduce  suitable  values  for  the  k-constants  that  were  omitted  from 
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(55)  etc.  The  preferred  values  of  kg,  k^,  k^  and  k  are  unchanged  from 

section  3,  but  k  ,  which  we  now  define  such  that  a.  contains  the  term 

-2  -  3  -  -  2  - 
k^a  /p  ,  has  the  untruncated  form  h  +  Jfe  cos  2 5  .  (This  gives  an  a  which  is 

much  more  convenient  than,  but  only  slightly  different  from,  the  Kozai  semi¬ 
major  axis,  for  which  the  k  is  given  by  (263);  for  the  absolute  constant  a'  , 

.  .  .  ,a  11 
introduced  in  section  3,  which  is  incidentally  the  semi-major  axis  of  Brouwer  , 

the  untruncated  k  is  given  by  (265),  as  follows  from  (64).)  We  also  require 

3 

k^  ,  not  previously  introduced;  if  this  is  defined  such  that  "+1^"  follows 

"3f  S^"  in  the  first  square  brackets  of  (61),  then  k^  is  required  to  be 

identical  to  k  -  it  will  be  observed  that  k  and  k  do  not  exactly  cancel 
0)  w  M 

in  U  (they  would  if  U  were  qio  +  M  instead  of  co  +  M  ) ,  the  effective 

_2 

constant  in  Uj  being  ^e  k^/(I  +  q)  . 

On  this  basis,  elaboration  of  (257),  (258)  and  (259)  leads  to  the  formulae 


ia  |fq2  cos  2u  +  he2  [2  -  (r/p)  (5  -  cos  2v)jj  , 

|f  [sin  2u  +  4e  sin  (u  +  to)]]  +  2h  (e/q2)  [4(3  -  4e2)  sin  v 

-  e  sin  2v  -  4$eq2  V(v,  e)]| 


(260) 


(261) 


^■re  sin  2i  ^2  sin  (u  +  v)  -  3  sin  5  -  ife  V(v,  e)  cos  u|  , 


(262) 


where  V(v,  e)  = 


(v  -  M  -  2e  sin  v) 


*  sin  2v  +  0(e) 

by  (50). 

Formulae  (260)  and  (261)  are  respectively  compatible  with  the  expressions 

for  dr  /p  and  du  +  dll  cos  i  given  by  Merson  on  page  8  of  Ref  7,  if  allow- 
s  s  s 

ance  is  made  for  the  different  values  of  k  and  k  employed;  Merson  was 

3  6 

working  with  the  unmodified  Kozai  elements  that  are  obtained  by  taking 


k  =  $h  (2  +  3e2  +  q2)  +  »fe2  cos  25 
3 

and 

kg  =  1 3 fi  |  5  +  2q2/(I  +  q)j  +  3f  cos  25  . 


(263) 

(264) 


Formula  (262)  effectively  combines  (but  with  k^  here  unity  instead  of  zero) 
Merson' s  expressions  for  dig  and  dilg  sin  i  . 

In  using  (260)  to  (262),  the  first  step  is  the  generation  of  mean  position, 
as  usual,  and  this  raises  no  problem,  since  untruncated  ft  and  w  are  already 
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available  from  and  S  ,  as  given  by  (67)  and  (78).  If  we  wish,  we  can 

continue  to  use  h  as  the  rate  of  change  of  U  ,  as  opposed  to  M  ,  if  we  make 

-2  .  -  . 
an  0(e  )  correction  to  and  ,  but  it  is  now  perhaps  more  natural  to 

associate  n  with  M  .  Then  n  no  longer  includes  C  and  becomes  the  standard 

6  J  1 

mean  motion  used  by  both  Kozai  and  Brouwer  ;  this  is  actually  identical  with 
our  n' ,  and  the  associated  (untruncated)  is  given  by 

kn  =  ^  (2  +  3l2>  +  cos  25  »  (265) 

as  opposed  to  the  previously  ’preferred’  value  given  by  (109).  There  is  a  con¬ 
sequent  change  in  p ^  ,  which  is  given  by 

i,  =  -  h  (1  -  3e2)  , 

if  k  is  standard  (or  by  -hq  if  k  for  the  unmodified  Kozai  semi-major  axis  is 
a.  a 

being  used),  instead  of  (162).  The  changes  in  k  and  p.  have  an  inevitable 
^  ^  * 
effect  on  p^  ,  even  if  this  is  still  defined  so  as  to  make  zero;  the  new 

value  is  given  by 

w2  =  -  J4f  (20  -  1  If)  , 


instead  of  (185)  . 

9  ON  SECULAR  AND  LONG-PERIODIC  PERTURBATIONS 
9 . 1  General  remarks 

The  original  intention  was  to  limit  this  paper  strictly  to  e-independent 

2 

J2  perturbations,  and  it  has  been  seen  that  the  resulting  expressions  for  short- 

periodic  perturbations  are  very  simple.  However,  for  an  eccentricity  of  0.01, 

.  .  2 
say,  a  low-altitude  satellite  experiences  long-term  perturbations  of  order  e 

that,  within  about  a  day  (since  this  corresponds  to  an  angular  motion  of  about 

100  radians),  are  of  the  same  order  of  magnitude  as  the  short-periodic  perturba- 

2 

tions  of  order  J.  .  Hence  the  paper  would  not  be  complete  without  some  reference 
z  2 

to  these  perturbations.  Once  we  go  to  order  e  ,  there  is  not  much  difficulty 
in  giving  the  full  formulae,  valid  for  any  eccentricity. 

The  long-term  e-independent  perturbations  are  purely  secular.  When 
e-dependent  terms  are  considered,  however,  the  long-term  perturbations  contain 
long-periodic  components,  trigonometrically  related  to  the  argument  of  perigee 
(u>) ,  as  well  as  purely  secular  components.  The  nature  of  these  perturbations 
is  frequently  misunderstood,  so  some  general  remarks  are  offered  before  the 
results  for  the  elements  a,  e,  i,  ft,  u>  and  M  are  given.  (There  is  no  need 


to  go  to  non-singular  elements  5,  o  and  U  ,  since  e  '  factors  do  not  appear 
in  the  long-term  perturbations.) 

Suppose  a  term  Kn  cos  kw  (with  K  not  necessarily  the  ^-related 
quantity  used  in  this  paper)  occurs  as  a  component  of  the  rate  of  change  of  some 
element,  C  say.  If  io  is  constant  (as  it  is  for  orbits  at  the  critical 
inclination)  we  can  integrate  to  get  a  secular  perturbation  Knt  cos  kw  .  If  w 
is  not  constant,  but  is  itself  essentially  secular  with  rate  w  ,  the  perturba¬ 
tion  can  be  written  in  the  form 

[sin  k£j  1  , 

*P  k5 

where  [....J  designates  "variation  over  the  period  of  integration",  te  it  is 
a  definite  integral.  The  two  results  are  compatible,  since  [sin  ku  J/ku3  tends 
to  t  cos  ku  ,  over  the  interval  t  ,  as  u  tends  to  zero.  If,  however,  6c 
is  written  with  an  arbitrary  lower  limit  of  integration  as  with  short-periodic 
perturbations,  in  particular  in  the  form  sin  ku/ku  ,  the  perturbation  has  an 
apparent  singularity  for  Si  zero. 

As  was  made  clear  by  Merson7,  trouble  is  avoided  if  we  keep  to  definite 
integrals.  Since  we  are  concerned  with  k  *  2  ,  we  denote  by  I£  and  Ig  the 
time  integrals  of  cos  2u  and  sin  2u  -  in  Ref  7  the  negatives  of  these  quanti¬ 
ties  were  denoted  C 2  and  S2  respectively. 

The  difficulty  is  largely  a  matter  of  terminology*.  The  terms  'first-order', 

'second-order',  etc,  are  unambiguous  when  used  to  describe  short-periodic  and 

2  .  ... 

secular  perturbations,  a  second-order  (J„)  perturbation  having  the  intuitive 

.  .  1  2 
meaning  that  the  perturbation  would  have  been  X  times  larger  if  J2  had  been 

X  times  larger.  It  is  not  so  simple  for  a  long-periodic  perturbation,  however, 

since  the  frequency  of  such  a  perturbation  increases  by  a  factor  of  X  when  its 

amplitude  increases,  thereby  reducing  the  increase  in  magnitude  of  the  long-term 

variation. 

It  seems  most  natural  to  take  the  order  of  a  long-periodic  perturbation 
from  the  power  of  J2  appearing  in  the  appropriate  rate-of-change  component,  ie 
in  C  for  element  C  .  With  this  philosophy,  there  are  no  long-periodic 

*  The  long-standing  division  of  opinion  over  the  problem  of  the  critical  inclin¬ 
ation  is  a  good  example  of  the  semantic  confusion  related  to  the  meaning  of 
'order' .  As  Allan  has  pointed  out^,  the  behaviour  of  certain  of  the  orbital 
elements  over  very  long  periods  of  time  is  indeed  a  libration  of  amplitude 
proportional  to  J2  ,  but  this  does  not  imply  a  singularity  -  the  effect  arises 
from  rates  of  change  that  are  still  only  of  order  j2  . 


100 


59 


perturbations  of  first-order  (though  J 2  is  exceptional  here  and  there  are 

such  perturbations  for  eta).  Second-order  perturbations  arise,  and  may  be 

-2 J  -2 

expressed  in  terms  of  K  I  and  K  I  .  Since 

c  s 

I£  =  [sin  2£5  J  /2S  , 

•  —  — 

where*  5  =  £Kn  (4  -  5f)  , 


K  I  behaves  as  a  first-order  perturbation,  and  is  regularly  so  described,  but 
in  the  present  paper  it  will  be  regarded  as  still  of  second  order.  Such  pertur¬ 
bations  are,  after  all,  not  merely  less  important  than  first-order  secular 
perturbations  -  they  are  also  less  important  than  second-order  secular 
perturbations.  (Their  long-term  magnitude  is  admittedly  much  greater  than  for 
short- periodic  perturbations  of  the  second  order,  but  this  is  no  more  than  the 
usual  more-important-if-longer-period  effect.) 


The  above  is  not  the  full  story,  however,  as  we  shall  encounter  some 
apparently  first-order  long-periodic  perturbations  that  are  really  of  the  third 

JL  ±  —  _  — 

order.  The  dominant  secular  perturbations  Q  and  w  ,  given  by  -  Kn  cos  i 
and  iKn  (4  -  5f)  respectively,  are  functions  of  e  (an  argument  of  K  )  and 
i  ,  so  that  they  have  induced  long-periodic  variation  -  thus  cos  i  is  equal  to 
an  initial  (hence  constant)  value,  cos  i^  ,  plus  a  term  given  by  -  6i^  sin  i^  . 
The  integral  of  K  <5i  is  itself  long-periodic  and  of  third  order,  but  because 

x,p 

two  integrations  are  now  involved  it  is  apparently  first  order,  and  particularly 


so  if  written  with  a  factor  E 


•-2 


Since,  as  we  shall  see,  6e„  and  <5i„  are 
*  ’  Up  Jlp 


multiples  of  I  ,  it  is  natural  to  introduce  the  notation  II  for  the  integral 
s  s 

of  this  quantity,  ie  for  the  double  time  integral  of  sin  2ui  . 


Two  further  preparatory  points  must  be  made.  First,  the  form  each  long- 
periodic  perturbation  takes  is  dependent  on  the  expression  for  the  first-order 
short-periodic  perturbation,  in  that  the  arbitrary  k-constant  may  be  chosen,  in 
different  ways,  as  a  function  of  u  .  Second,  there  is  an  entirely  different  way 
of  looking  at  second-order  long-term  variation,  and  this  was  the  approach  adopted 
in  Refs  2  and  4,  in  particular.  It  involves  the  derivation  of  each  element's 
variation  over  a  complete  revolution  of  the  satellite,  to  eliminate  the  short- 
periodic  effect,  but  now  the  interpretation  of  "complete  revolution"  is  crucial  - 


is  just  the  first-order  secular  rate,  ie  Kw 


Here 

by  the  second-order 


K2  E2 


but  the  time  derivative  of 


We  can  correct  this 
C  is  no  longer 


actually  constant,  now  that  we  consider  w  to  have  long-periodic  as  well  as 


secular  variation.  The  quantities 
long-periodic  variation. 


Ic,  Is  and  IIS  necessarily  neglect  the 


the  variation  from  perigee  to  perigee,  for  example,  is  different  from  the  varia¬ 
tion  from  node  to  node,  not  because  perigee  and  node  are  different  points  but 
simply  because  there  is  a  slippage  (due  to  55  ) ,  as  between  successive  perigees 
and  successive  nodes. 

9.2  Perturbation  in  a 

There  is  no  secular  perturbation  in  the  semi-major  axis,  nor  is  there  any 
second-order  long-periodic  variation  so  long  as  the  w-dependent  part  of  the 
'constant'  k^  is  specifically  |fe  cos  2w  .  This  is  the  case  not  only  for  the 
three  ka's  that  have  been  quoted  -  the  preferred  value  of  this  paper,  the 
Brouwer  value  associated  with  a'  and  the  Kozai  value  -  but  also,  it  appears, 

•  •  pc 

for  the  k  implied  by  two  of  the  papers  (the  one  by  Berger  and  Walch  and 
a  26  " 
the  one  by  Kinoshita  )  considered  in  section  5.  (The  value  of  this  k  is 

-  3  a 

less  than  the  Kozai  value  by  hq  .)  Given  an  arbitrary  k  ,  however,  we  have 

cl 

a  long-periodic  perturbation  given  by 

<5a£p  =  Kaq-2  |[2f£2  cos  2S  -  k_J  .  (266) 

As  remarked  in  section  9.1,  there  is  an  apparent  long-periodic  effect  if 
we  only  consider  the  value  of  a  at,  say,  perigees  or  ascending  nodes.  For 
perigees  the  effect  is  given  by 

=  -  K2naq“2  (1  +  e)3  f  (4  -  5f)  Ig  , 

as  indicated  by  equation  (181)  of  Ref  4. 

9 . 3  Perturbation  in  e 

There  is  again  no  secular  perturbation,  but  there  is  a  long-periodic 
perturbation  given  by 

6es.p  =  Hi9  f  (I4  "  15*>  !s  +  i  Ke  [3f  cos  2to  -  kfi  1  ,  (267) 

where  the  second  term  vanishes  for  both  our  standard  kg  and  the  value  implied 
by  use  of  Kozai  elements.  However,  Refs  25  and  26  involve  a  value  of  k  that 

l --2  _  _  e 

is  less  than  Kozai's  value  by  l$fq  (1  +  2q)(l  +  q)  cos  2S  ,  whence  it  follows 
that,  to  0(i)  ,  their  expression  for  5e  should  involve  a  factor  2(13  -  15f) 
in  place  of  the  factor  (14  -  15f)  in  the  standard  expression.  This  is  con¬ 
firmed  by  BW  112.15  ,  in  the  notation  of  section  5.2,  referring  to  Ref  25. 
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9.4  Perturbation  in  i 

2 

A  simple  relation  exists  between  perturbations  in  a,  e,  i,  since  p  cos  i 
is  constant,  as  remarked  in  section  5.3.  This  constancy  applies  to  any  potential 
field  that  is  symmetric  about  the  z-axis,  ie  to  the  potential  associated  with 
any  zonal  harmonic.  Assuming  6a.  to  be  zero,  it  follows  that 

6is.p  =  cot  1  •  (268) 

In  fact 

6i£p  =  5sK2ne2  sin  2i  (14  -  15f)  Ig  -  JK  sin  2i  [k.]  ,  (269) 


where  the  second  term  vanishes  for  our  standard  k^(  =  1)  and  the  Kozai 

k. (  =  0)  .  However,  Refs  25  and  26  involve  a  value  of  k.  that  exceeds  Kozai' s 

value  by  ie  (1  +  2q)(l  +  q)  cos  2m  ,  whence  it  again  follows,  to  0(e)  ,  that 

their  expression  for  6i.  should  involve  a  factor  2(13  -  15f)  in  place  of 

_  *P 

the  factor  (14  -  15f)  in  the  standard  expression.  This  is  confirmed  by 
BW  112.18. 


9.5  Perturbation  in  Q 

The  complete  second-order  long-term  perturbation  may  be  obtained  from  the 
following  expression,  based  on  Ref  7,  for  the  rate  of  change  of  the  mean  element: 


||  =  -  Kn  cos  I  jl  +  i(“/n)  +  J^Kq"2  [4(15  -  19f) 

+  e2  (4  -  9f)  -  e4  (4  +  5f)  +  e2  cos  25  ((14  -  3f) 

-  2e2(7  -  15fj)  -  48  ka  +  12S2  kfi  -  12fq2  k£  -  36kn]^.  (270) 

We  consider  the  secular  component  first.  With  our  standard  values  of 

k  ,  k  and  k.  ,  and  k  given  by  (265),  we  get  (of  (137)  if  k  is  from  (109)) 
sc  i  n  n 

fl2  =  jVn  cos  I  {4  (3  -  5f)  -  e2  (4  +  5f)|  .  (271) 

With  the  Kozai  k' s,  on  the  other  hand,  we  get 

n2  =  jVn  cos  i  |48hq  -  4(9  -  lOf)  -  e2  (4  +  5f)}  , 

as  given  (effectively)  in  Refs  6  and  7.  Finally,  with  the  k's  of  Refs  25  and 

26,  we  get 


& 

1 


t 


3 

D 


fi2  -  -  jVH  cosi  |24hq  +  4(9  -  1  Of )  +  e2  (4  +  5f)| 
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in  agreement  with  BW  110.20-24  and  with  the  sixth  formula  of  KB-53  (notation 

of  section  5.4,  in  reference  to  Ref  26)  -  since  Ref  25  uses  (R/a)^  rather  than 
4 

(R/p)  as  (effectively)  a  factor  of  K  ,  the  bracketed  expression  here  must  be 
——8 

multiplied  by  q  before  agreement  is  evident. 


Turning  to  the  long-periodic  component,  we  can  write  the  basic  second-order 
perturbation  as 


ftp(b) 


-  j^K^ne^  (7  -  15f)  cos  i  1^  -  JK  cos  i  fk^J  +  4"K2nq  2  cos  i  x 
4  [4  (ka)^  ~  3fe2  cos  2JlQ  -  4e2  [(kfi)^  -  3f  cos  2w3 

+  4fq2  (k.)w+  3[4<kn>w  ~  3fe2  cos  2Ci]l  dt  ,  (272) 


where  (k  )  denotes  the  iv-dependent  (long-periodic)  part  of  k  ,  etc.  But 

3.  Cl)  ** 

in  addition  we  now  have  an  induced  (third-order)  perturbation,  as  explained  in 
section  9.1.  On  the  assumption  that  the  first-order  secular  perturbation  is 
computed  from 


=  -  Kq  Eq  cos  ig 


where  zero-suf f ices  have  been  added  to  indicate  that  we  use  initial  values  of 
the  mean  elements,  our  additional  perturbation  is  given  by 


^^Hp(ind) 

from  this,  using  (268), 


-  Kn  cos  i 


dt 


we  may  write 


^^Ap(ind) 


-  5Kneq  2 


cos  1 


i/«e 


£p 


dt  . 


Hence  our  standard  expression  for  the  combination  of  ^^  ^)  an<* 


6Vind)  “  giV6n  by 


£p 


=  -  5*3 K2ne2  cos  1  ^2(7  -  I5f)  I£  -  5K5f  (14  -  15F)  IIg 


(273) 


Use  of  the  k's  of  Refs  25  and  26,  on  the  other  hand,  including 
^e  (1  +  2q) (1  +  q)  sin  22  for  k^  ,  leads  to 

6n«,p  =  "  ^K2ne2  cos  I  |[2(I3  -  30f)  +  0(e2)]  lc  -  lOicHf  (13  -  15f)  II  ^  ; 


if  (by  taking  indefinite  integrals  instead  of  evaluating  Ic  and  II s  -  see  the 
discussion  in  section  9.1),  we  interpret  this  as  a  first-order  perturbation, 
we  get 
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K  e2  cos  I  (52  -  120f  +  75f2) 
12  (4  -  5f)2 


-4 

sm  2u  +  0(e  ) 


which  tallies  with  BW  112.20  . 

9.6  Perturbation  in  10 

The  complete  second-order  long-term  perturbations  may  be  obtained  from  the 
following  expression,  based  on  Ref  7,  for  the  rate  of  change  of  the  mean  element: 


=  iK5  {(4  -  5f  )  -  H^/H)  Sk  /3(o  +  ^Kq“2  [2(288  -  676f  +  395f2) 
at  \  a)  *- 

-  e2 (40  +  4f  -  65f2)  -  I4 (56  -  36f  -  45f2)  -  2  cos  2u  (20(5  -  6f) 


-  e2 (28  -  4f  + 


1 5f 2)  +  e4  (28  -  158f  +  135f2))-  24(4  -  5f)  x 


x  (4k,  -  e2k,  +  3k  J  -  240f(l  -  f)q2k.]J. 


n 


(274) 


We  consider  the  secular  component  first.  With  our  standard  values  of 

k  ,  k  and  k.  ,  and  k  given  by  (265),  we  get  ( cf  (216)  if  k  is  from  (109)) 
<*  e  1  n  n 

•>}  .  (275) 

With  the  Kozai  k's  ,  on  the  other  hand,  we  get 

u2  =  -  tfeH  |96hq  (4  -  5f)  -  2  (192  -  412f  +  215f)  -  £2  (56  -  36f  -  45f2)}  , 


h2  *  -  jj1? n  ^2f  (4  +  25f )  -  §2  (56  -  36f  -  45f2' 


as  given  (effectively)  in  Refs  6  and  7.  Finally,  with  the  k's  of  Refs  25  and 
26,  we  get 

32  =  9*sn  |48hq  (4  -  5f)  +  2(192  -  4 1 2f  +  215f2)  +  e2  (56  -  36f  -  45f2)}  , 

in  agreement  with  BW  110.25-29  and  with  the  fifth  formula  of  KB-53;  as  with 

•  _g 

*  the  bracketed  expression  must  be  multiplied  by  q  to  make  agreement  with 
Ref  25  evident. 


Turning  to  the  long-periodic  component,  we  can  write  the  basic  second-order 
perturbation  as 

6%(b)  “  “  *VK2n  {2f  (14  -  15f)  -  e2  (28  -  158f  +  135f2)|  Ic 

-+K  Ik^  +  3f  sin  u|  -  t^K2^-2  J j(4  -  5f)  x 
x  (4Q4  (ka)uj  -  3fl2  cos  23]  -  4I2  [(k  )  -  3f  cos  25] 

+  3 [4 (k  )  -  3fe2  cos  23])  +  40f(l  -  f)q2(k.)  I  dt  ,  (276) 

n  GO  ■y  1  Ok)  f 
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using  the  same  notation  as  for  ^£p(fc)  '  a^so  ^ave  an  induced  perturbation 
that  arises  from  computation  of  the  first-order  secular  perturbation  using 


S1  =  *K0  H0  (4  "  5f0)  ‘ 


This  additional  perturbation  is  given  by 


n  2eq 


5u>„  , .  ...  =  Kn 

£p(ind) 


this  gives,  using  (268), 

5u)„ 


(4  -  5f)  /fie  dt  -  5  sin  i  cos  i 

x,p 


/%dt] ; 


=  Kiieq'2  (13  -  15f) J Se„  dt 


£p(ind)  - -  --'J—t  p 

Hence  our  standard  expression  for  the  combination  of  6w 
is  given  by 

;2- 


„  ...  and  5u>„  .. 
Hp(b)  £p(ind) 


fcp 


48  I 


Zn  j[2f  (14  -  1 5f )  -  e2  (28  -  158f  +  135f2]  Ic 
+  2Kne2f  (13  -  15f)(14  -  I5f)  II 


}  • 


(277) 


Use  of  the  k’ s  of  Refs  25  and  26,  on  the  other  hand,  including 


$[3f  +  4(1  +  2q)  ( 1  +  q)~2 (f  +  e2  -  2e2f)]  sin  2<H 


for  k  ,  leads  to 

(D 


fill) . 


=  ”  {C2f  <13  ”  15f>  "  < 2 6  -  1 55f  +  140f *■) 

+  0(e4)]  I  +  2Kne2f  (13  -  15f)2  II  1  . 

J  c  s ) 

If  we  take  the  apparent  first-order  perturbation  that  results,  we  get 

K  sin  2u>  (2~  (52  _  ]25j  +  75f 2)  _  g2  (104  _  412f  +  555f2  -  250f3)  +  0(e4)|-  ; 


24  (4  -  5f)2 


-  -  4  2-2 

this  conforms  with  BW  112.23-24  if  we  remember  to  replace  K  by  Kq  (1  -  e  ) 

-2  -  -2  _•» 
so  that  the  coefficient  of  e  changes  to  -  (104  -  620f  +  I055f  -  550f  )  . 

9.7  Perturbation  in  M 


Merson  has  shown  that  the  rate  of  change  of  M  is  given  by 

-  n’  -  j4K2nq3  {(8  -  8f  -  5f2)  -  2f  (14  -  15f)  cos  23  j 

+  i Kgj  (31^/30)  +  6f  cos  23)  ,  (278) 
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where  the  final  term  has  been  added  here  to  cover  the  possibility  of  a  non¬ 
standard  value  of  k^  .  This  expression  is  'standard'  from  the  point  of  view  of 
the  present  paper. 

To  generalize  from  n'  ,  which  is  the  absolute  constant  introduced  in  1 

section  3,  to  an  arbitrary  n  would  require  a  deeper  analysis  than  will  be 

embarked  upon  here,  since  it  would  involve  the  derivation  of  the  values  of  k„ 

2n 

(see  section  4.8)  appropriate  to  n'  and  the  selected  ii  .  However,  the  diffi¬ 
culty  is  not  too  great  in  going  from  n'  to  the  n  of  Ref  25,  so  long  as  it  is 
taken  only  to  0(e).  The  reward  for  this  is  an  overall  check  between  Refs  7  and 

25  and  the  formula  for  k„  in  section  4.1 

2a 

The  values  of  k^n  ,  appropriate  to  n'  and  n  (the  n  of  Ref  25  will 

be  understood  in  what  follows),  can  be  derived  from  the  values  of  k„  ,  with  p, 

L  a  1 

and  u  (see  section  4.8)  taken  as  zero  in  both  cases.  Thus  (161)  gives  k  =  k 
z  n  3 

and  then  (168)  gives 

k2„  '  *  (5f2  *  ,0ka  '  8k2a>  ' 

2  3  ”2 

For  n'  (given  by  n'  a'  =  y  )  we  have  (untruncated)  k  =  $h(l  +  lje  )  + 

—  2  _  a  _ 

+  i fe  cos  2C  ,  as  quoted  in  section  8,  and  k^a  given,  to  0(e)  ,  by  (117). 

This  leads  to  what  we  may  designate  k_  (n')  ,  an  expression  for  which  has 

zn  2  3 

already  been  given  -  see  (188).  For  n  (which  satisfies  n  a  =  y)  ,  on  the 

2 

other  hand,  both  k  and  k  are  0(e  )  so  that,  to  0(e)  , 

SL  Z  ci 

k2n(n)  =  Hf2  . 

But  the  k„  's  are  defined  such  that 
Zn 

n'  |l  -  ljK'  kn(n')  +  K'2  k2n(n')} 

and 

n  |l  -  1|K  kn(n)  +  K2  k2n(n)| 

are  identical  to  0(e)  ,  since  they  represent  the  same  mean  value  of  the  osculating 

mean  motion,  n  .  In  making  the  identification  we  can  take  k^  to  be  zero,  in 

2 

(188),  since  it  is  0(1  )  in  Ref  25.  Then  we  obtain,  to  0(e)  , 

n'  =  n  +  K’  n'  h  +  j^K2  n  (40  -  48f  +  30f2)  . 

But  K'  n'  must  be  replaced  by  K  n  in  the  first-order  term  of  this  relation, 
where  in  fact 

-A 
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K'  n'  =  Kn  (1  +  2^-ich)  +  0(K2  e2)  , 

since,  to  0(e)  ,  K'  =  K(  1  +  l-^Kh)  and  n'  -  ii(l  +  Kh)  .  Hence  we  get 

n'  =  n  +  Knh  +  ?*gK2  n  (152  -  384f  +  282f2)  +  0(e2)  .  (279) 

_2 

On  combining  the  K  terms  in  (278)  and  (279)  ,  we  see  that ,  residual  to  the  n 
of  Ref  25  and  the  first-order  term  Kiih  ,  dM/dt  contains  the  secular  component 

Aic2  n  (144  -  376f  +  287f2)  . 

This  tallies  with  BW  110.30  (Ref  25),  but  not  with  the  fourth  formula  of  KB-53 

(Ref  26).  However,  the  apparent  error  in  second-order  M  ,  as  given  by  Ref  26, 

has  already  been  remarked  (sections  5.2  and  5.4)  and  is  resolved  by  allowance  for 

an  effective  non-zero  k„ 

2a 

The  long-periodic  perturbation  is  easy  to  express  in  the  appropriate  form 
for  checking  against  Ref  25.  The  appropriate  k^  is  -  ■J’f  |(3  +  2(2  +  e2)  x 
(1  +  2q) / ( 1  +  q2)|  sin  23  ,  whereupon  the  last-term  of  (278)  yields 

nf  (4  -  5f)^6  -  7e2  +  0(e4))  cos  2m  . 

This  leads  to  a  perturbation  of  the  form 

^K2nf  ^2(13  -  15f)  -  5e2  (7  -  8f)  +  0(e4)^  Ic  ;  (280) 

this  conforms  with  BW  112.27  and  BW  112.28  if  we  remember  to  replace  K  by 
-4-2-2  2 

Kq  (1  -  e  )  ,  so  that  the  coefficient  of  §  changes  to  (17  -  20f)  ,  and 

also  note  that  BW  112.27,  28,  29  and  30  all  contain  an  unnecessary  factor, 

(4  -  5f)  ,  which  may  be  divided  out. 

10  DISCUSSION  AND  CONCLUSIONS 

The  main  goal  of  the  paper  has  been  the  derivation  of  the  principal  ( ie 
.  2 

e-independent)  ^  perturbations  of  a  satellite  orbit  by  an  elementary  approach, 
and  the  expression  of  the  results  in  as  compact  a  set  of  formulae  as  possible. 

The  elementary  approach,  assuming  nothing  beyond  the  planetary  equations  of 
Lagrange,  has  inevitably  involved  some  long  and  tedious  algebra,  but  simple 
formulae  were  eventually  obtained  by  expressing  short-periodic  perturbations  in 
terms  of  a  rotating  system  of  cylindrical  polar  coordinates. 

The  introduction  of  the  paper  referred  to  the  sophistication  of  the  litera¬ 
ture  and  to  the  difficulty  in  picking  out  the  main  results  from  a  mass  of 


mathematics.  It  is  therefore  necessary,  to  avoid  the  same  charge,  that  the 

present  author  summarize  his  main  results,  either  directly  or  by  reference  to 

the  equation  numbers.  This  will  first  be  done  on  the  basis  of  the  assumption 

made  during  most  of  the  analysis,  namely,  that  eccentricity  is  of  the  same 

2  2 

order  of  magnitude  as  Jj  ,  so  that  J^e  ,  as  well  as  J^e  ,  can  be  neglected; 

there  is  also,  of  course,  the  tacit  understanding  that  the  results  are  not  to  be 

3 

used  over  such  a  long  period  of  time  as  would  make  J  perturbations  intolerable, 

2  2  1 

as  well  as  the  long-term  and  perturbations.  Secondly,  complete 

first-order  results  are  indicated,  permitting  considerable  relaxation  on  the 

smallness  of  eccentricity.  Finally,  to  provide  results  that  are  valid  over  an 

2 

extended  period,  the  complete  J ^  formulae  for  long-term  development  are  indi¬ 
cated,  involving  the  special  treatment  of  long-periodic  perturbations. 

The  starting  point  is  assumed  to  be  a  set  of  mean  elements  at  epoch,  viz 
V  eQ,  iQ,  fiQ,  and  MQ  ,  whence  lQ,  nQ  and  UQ  by  (1),  (2)  and  (3); 

values  of  y,  and  R  are  also  assumed,  whence  also  K  following  (21). 

The  actual  generation  of  these  mean  elements,  given  osculating  elements  (or  the 
coordinates  of  position  and  velocity),  would  not  be  entirely  trivial,  but  if  no 
other  recourse  were  available  they  could  be  adjusted  iteratively  so  that  the 
perturbations  led  back  to  the  osculating  elements.  The  significance  of  'mean 
elements'  was  discussed  in  section  2.4,  and  the  paper  contains  many  general 
formulae  which  can  be  evaluated  in  terms  of  particular  definitions  of  these 
elements . 

For  our  first,  and  most  restricted,  derivation  of  satellite  position  at 
time  t  ,  we  start  with  a,  e  and  i  ,  unchanged  from  a^,  5^  and  iQ  ,  and 
hence  derive  n  from  the  version  of  Kepler's  third  law  expressed  by  (256). 

The  values  of  ft,  u  and  U  are  given  by 


n 

=  -  Kn  cos  i 

{' 

-  |K  (3  -  4f ) |  t 

,  (281) 

(0 

g 

+ 

o 

13 

II 

(4  - 

5f) 

t 

(282) 

and 

U 

=  °o  +  Ht 

» 

(283) 

where  the  expressions 

•  • 

for 

and 

“l 

,  given  by  (67), 

(137)  and  (78),  have 

been  employed. 

The  mean  elements  lead  at  once,  by  the  standard  algorithm  of  section  2.5, 
to  r  and  u  .  Then  (88),  (102),  (100),  (228),  (246)  and  (252)  can  be  combined 
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r  =  r  +  &Kaf  |cos  2u  -  nK  [f  cos  4u  +  2(26  -  31  f)  cos  2u]]|  ,  (284) 

u'  =  u  +  iS'K  'jf  sin  2u  +  4e  [6h  sin  v  +  f  sin  (u  +  to)3 

-  iKf  [f  sin  40  -  ( !  9  -  20f)  sin  2u]|  (285) 

and 

c  =  ^Kr  sin  2i  je  [2  sin  (u  +  v)  -  3  sin  uQ  -  iKf  sin  3u|  ,  (286) 

to  give  the  special  cylindrical  coordinates.  Finally,  x,  y  and  z  are  given 
by  (45),  r'  being  identified  with  r  . 

Expressions  for  osculating  elements,  if  these  are  desired,  are  necessarily 
more  complicated  than  (284)  to  (286).  The  element  £  is  given  by 

C  =  l  +  SC  , 

where  first-order  contributions  to  <5£  ,  for  £  equal  to  a,  i,  12,  £,  q  and  U  , 
are  given  by  (65) ,  (66)  (68) ,  (70)  (72)  and  (74) ,  and  second-order  contributions 
by  (122),  (128),  (136),  (142),  (147)  and  (180);  the  latter  involve  preferred 
values  of  the  k-constants  employed  throughout  this  paper,  the  generalized  expres¬ 
sions  being  given  by  (121),  (127),  (134),  ( 1 4 1 ) ,  (146)  and  ( 179 )  -  correspondingly 
•  • 

generalized  0^  anc*  U 2  are  8i-ven  by  (133)  and  (181). 

If  the  first-order  effects  are  to  be  covered  completely,  ie  without 
e-truncation,  it  is  perhaps  more  natural  to  define  n  so  that  (283)  is  replaced 
by 

M  =  MQ  +  5t  ,  (287) 

though  the  continued  absence  of  a  second-order  component  of  C  means  that  the 
second-order  part  of  n  still  really  relates  to  U  .  This  implies  a  change  of 
,  with  a  new  version  of  Kepler's  third  law,  given  by 

n2  a3  =  u  |l  -  Kh  (1  -  3e2)  -  Aric2  f  (20  -  llf)|>  (288) 

instead  of  by  (256),  where  the  second-order  component  now  comes  from  the  last 

equation  of  section  8.  The  change  in  k  affects  (281),  which  has  to  be 

n 

replaced  by 

a  =  Qq  -  Kn  cos  i  {l  -  £k  (3  -  5f)[  t  ,  (289) 

and  the  first-order  components  of  r,  u'  and  c  are  now  as  given  by  (260)  to 
(262)  and  net  as  in  (284)  to  (286), 
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If  the  long-term  second-order  effects  are  to  be  represented  completely, 
with  an  explicit  second-order  component  of  S3  in  particular,  then  to  identify 
n  with  M  we  require  yet  a  further  version  of  Kepler's  third  law,  viz 


-2  _3 

n  a 


"  {'  -  Wo  (l  -35o>  -**oFo  «  -9!o>}  • 


(290) 


in  which  zero-suff ices  appear  because  e  and  i  are  no  longer  constant.  The 
expressions  for  the  long-term  variation  of  the  mean  elements  are  as  follows: 


a  =  a. 


(291) 


e  = 

-  *K^n 

■o$o  <14  -  '“o>  ’s  • 

(292) 

1  =  To 
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+  rtKQn 

*0  si”  210  04  '  15 V  Is 

9 

(293) 

n  =  0^ 

-  K0n  cos  I0  1  -*Kq[4(3  -  5f0) 

-  4(i  *  “0>]  | 4 

-  ^K^n 

*0  cos  T0  {2(7  ‘  l5?0nc 

-  V,"4  -  ’“o’11.}  • 

(294) 

O 

1  3 

II 

1  3 
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]|* 

-  ifa  K^n 
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(296) 


whe'e  I  ,  I  and  II  are  the  definite  integrals  of  cos  2 (ui_  +  wt)  , 
css  0 

sin  2  (u)q  +  Sgt)  and  Ig  ,  respectively. 

-2-2 

No  kq6q  contribution  is  required  in  (290)  because  we  have  not  taken  the 

short-periodic  perturbations,  in  a  in  particular,  to  O(K^e)  ,  let  alone 
2  2 

0(K  e  )  ;  hence  the  second-order  constant  k„  ,  introduced  in  section  4.2,  can 

z  a 

be  regarded  as  including  whatever  component  is  necessary  to  eliminate  any 
O(K^e^)  term. 
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